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ABSTRACT 

Acoustic  fatigue  is  the  result  of  high  frequency  lateral  vibration  of  an  ai  rcraft  panel  due  to  time 
varying  pressurewaves  caused  by  engineand/  or  aerodynamic  effects.  For  example,  acoustically 
induced  cracks  have  been  recorded  in  the  lower  external  su  rface  of  the  nacelle  skin  of  the  F/  A-18 
aircraft  and  aftfuselage.  In  the  case  of  the  inlet  nacelle  overall  sound  pressurelevelsoftheorder 
of  172  dB  havebeen  recorded.  Attempts  to  repair  these  cracks  by  applying  standard  methods  of 
bonded  repai  r  dev  el  oped  for  i  n-pl  ane  I  oads  were  made.  FI  owever,  the  cracks  conti  nued  to  grow 
at  a  si  mi  I  ar  rate.  W  h  i  I  e  th  e  repai  r  of  cracked  ai  rcraft  str  u  ctu  res  su  bj  ected  toin-planeloadsusing 
bonded  repairs  has  resulted  in  con  si  derable  aircraft  I  ifeti  me  extension  and  hence  cost  savings,  the 
use  of  bonded  patches  to  repair  panelswith  acoustically  induced  cracks  (acoustic  fatigue)  isonly 
recent.  I  n  this  report  a  generic  design  procedure  is  presented  for  the  repai  r  of  panels  contai  ni  ng 
acoustically  induced  cracks  using  constrained  layer  damping  (CLD).Theapplication  of  bonded 
repairs  to  acoustically-induced  cracks  requires  analytical  tools  that  take  into  account  high 
frequency  out-of-planevibrati  on. Theanalyti  cal  tools  described  in  the  report  will  enablethe  rapid 
desi  gn  of  effecti  ve  repai  rs  usi  ng  cl  osed  form  sol  uti  ons.  A  potenti  al  outcome  for  thi  s  study  coul  d 
be  application  to  the  repair  of  the  batch  11  F/  A-18  aircraft  which  suffer  from  acoustic  fatigue 
cracking  on  the  aftfuselage.  Over  a  period  often  years  this  could  result  in  cost  savings  of  one 
order  of  magnitude  i  improvement  over  those  for  metal  I  i  c  repai  rs. 
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Generic  Design  Procedures  for  the  Repair  of 
Acoustically  Damaged  Panels 


Executi  ve  Su  mmary 

C  u  rrentl  y  a  batch  of  el  even  RA  A  F  F/  A  - 18  ai  rcraft  are  d  evel  opi  ng  acousti  c  fati  gue  cracks  i  n 
the  rear  fusel  age  area  as  a  resu  1 1  of  hi  gh  frequency  I  ateral  vi  brati  ons  d  ue  to  ti  me  vary i  ng 
pressure  waves  caused  by  engine  and/  or  aerodynamic  effects.  The  repair  of  cracked 
aircraft  structures  subjected  to  in-plane  loads  using  bonded  repairs  has  resulted  in 
considerable  aircraft  life  time  extension  and  hence  cost  savings.  Flowever,  the  use  of 
bonded  patches  to  repair  panels  with  acoustically  induced  cracks  (acoustic  fatigue)  is  a 
recent  dev  el  opment.  I  n  thi  s  report,  a  generi  c  desi  gn  procedure  i  s  descri  bed  for  the  repai  r 
of  panels  containing  acoustically  induced  cracks.  This  procedure  consists  of  a  simplified 
dynamic  analysis  together  with  crack  bridging  theory.  Furthermore,  the  solution  is  also 
appl  i  cabl  efor  plates  contai  ni  ng  cracks  on  the  boundary.  I  n  additi  on  the  procedure  makes 
use  of  a  closed  form  solution  to  determine  the  effect  of  residual  thermal  stresses  and 
softening  of  the  adhesive.  It  has  been  shown  that  the  analysis  for  crack  bridging  can 
predict  the  behaviour  of  plates  simply  supported  on  all  edges  or  clamped  edges. 

The  generic  design  proceduremodelstheviscoelasticdamping  material  and  boron/  epoxy 
constraining  layer  in  order  to  damp  out  vibration  of  the  panels.  As  a  result  this  generic 
design  procedure  together  with  damping  data  will  aid  in  the  design  of  highly  damped 
repairs  to  acoustically-induced  cracked  panels.  The  generic  design  procedure  has  been 
val  i  d ated  usi  ng  f  i  n i  te  el  ement  anal  ysi  s.  The  goal  of  th i  s  w ork  i  s  eventu al  I  y  to  i  ncorporate 
the  design  method  into  the  RAAF  Engineering  Standard  DEF(AUST)9005-A,  to  enable 
rapid  design  of  such  repairs. 

It  has  been  found  that  the  use  of  a  highly  damped  repair  is  a  low  risk,  low  cost  repair  that 
could  be  used  on  the  aft  fuselage  of  the  F/  A-18  to  extend  the  life  of  this  component 
beyond  6000  hours.  Furthermoreover  a  period  of  10  years,  the  cost  savings  are  an  order  of 
magnitude  improvement  over  those  for  metallic  repairs. 
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1.  Introduction 

The  repai  r  of  cracked  ai  rcraft  structures  subjected  to  i  n-pl ane  I oads  usi  ng  bonded  repai  rs  has 
resulted  in  considerable  aircraft  lifetime  extension  and  hence  cost  savings  [1].  However,  the 
use  of  bonded  patches  to  repair  panels  with  acoustically  induced  cracks  (acousticfatigue)  isa 
recent  development.  In  this  report,  a  generic  design  procedure  for  the  repair  of  panels 
contai  ni  ng  acousti  cal  I y  i  nd uced  cracks  will  be  descri  bed .  A  cousti  c  fati  gue  i  s  a  resul  t  of  hi  gh 
frequency  lateral  vibration  in  aircraft  panels  due  to  time  varying  pressure  waves  caused  by 
engine  and/  or  aerodynamic  effects.  For  instance,  acoustically  induced  cracks  have  been 
recorded  i  n  the  I  ower  external  surface  of  the  nacel  I  e  ski  n  of  the  F/  A-18  ai  rcraft,  as  show  n  i  n 
Figure!  In  this  case,  overall  sound  pressure  levels  of  theorder  of  172  dB  have  been  reported, 
[2]  and  [3].  Attempts  were  made  to  repair  these  cracks  by  applying  standard  methods  of 
bonded  repair  [4]  developed  for  in-planeloads,  however,  thecracks  continued  to  grow.  It  was 
found  [5]  that  whi  I  ethe  boron  fibre  composite  bonded  repai  rsdid  reduce  the  stress  intensity 
at  the  crack-ti  p,  the  hi  gh  number  of  cyd  es  sti  1 1 1  ead  to  si  gnifi  cant  crack  growth  of  up  to  2  mm 
per  flight  hour.  The  application  of  bonded  repairs  to  acoustically-induced  cracks  requires 
analytical  toolsthat  can  takeinto  account  high  frequency  out-of-plane  vibration. 


Aft  fuselage  cracks 

Figurel:  Location  of  acoustic  fati  gue  cracking  in  F/A  -18  fusel  age  pan  els 

Si  nceit  is  known  thattheamplitudeof  vibration  is  inversely  proportional  to  the  square  root  of 
thedamping  ratio,  a  combi  nation  of  damping  and  stiffness  may  reduce  the  crack  growth  rate 
significantly.  The  application  of  highly  damped  repairs  (Durability  Patches)  to  acoustically 
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damaged  panels  has  been  proposed  by  Rogers  et  al  [6].  Also,  a  highly  damped  repair  was 
applied  [7]  to  the  vertical  fin  of  the  F- 15  aircraft,  as  part  of  a  test  program.  (More  recently,  a 
series  of  repairs  have  been  designed  and  tested  in  theAFRL  progressive  wave  tunnel  [8, 9]). 
Also  Boeing  [10]  havedeveloped  procedures  for  the  use  of  the  Dosimeter,  and  methodsfor 
optimal  design  of  damped  repairs.  Repairs  have  been  madeto  B52  fusel  age  skin  panels,  F-15 
access  panel  and  C 130 flap  panels,  and  theDosi meter  has  been  used  fortheassessementof  the 
effecti  veness  of  the  repai  r  [  10] .  The  D  osi  meter  bei  ng  a  smal  I  autonomous  devi  ce,  carri  ed  i  n  an 
aircraft,  and  designed  to  measure thethermo-acoustic  environment  of  an  aircraft  skin  plate. 

Inthisreport,  a  generic  design  procedure  (based  on  constrai  ned  layer  damping  (CLD))will  be 
presented .  Thi  s  w  i  1 1  i  nvol  ve  the  sel  ecti  on  of  a  dampi  ng  materi  al  and  the  geometry  of  the  C  L  D 
in  order  to  prevent  the  growth  of  acoustically-induced  cracks.  This  work  makes  use  of  the 
closed  form  solution  for  the  stress  intensity  factor  developed  by  DSTO  [11]. 

2.  Overview  of  design  procedure 

The  use  of  C  L  D  contai  ni  ng  a  vi  scoel  asti  c  materi  al  i  s  wel  I  know  n  and  i  s  an  effecti  ve  method  of 
di ssipating  the  energy  dueto  flexural  vibration.  As  shown  in  Figure2,  CLD  involves  a  plate, 
vi  scoel  asti  c  materi  al  and  a  constrai  ni  ng  I  ayer.  I  n  the  present  context,  the  constrai  ni  ng  I  ayer  i  s 
effectively  the  repair  patch.  The  bending  of  the  structure  involves  shear  deformation  of  the 
vi  scoel  asti  c  materi  al .  D  ampi  ng  occu  rs  si  nee  the  vi  scoel  asti  c  materi  al  i  s  a  pol  ymer  compri  sed 
of  long  chain  molecules  which  when  deformed  results  in  friction  of  these  molecules.  A 
measure  of  the  effectiveness  of  the  CLD  is  the  loss  factor,  77 .  For  a  vibrating  plate,  the  loss 
factor  may  be  defined  as  the  ratio  of  the  energy  dissipated  per  radian  to  the  total  strain 
energy.  The  viscoelastic  material  has  a  very  high  loss  factor  in  comparison  to  metallic 
components.  Also  the  constraining  layer  should  be  very  stiff  in  order  to  promote  shear 
deformati  on  of  the  vi  scoel  asti  c  materi  al ,  i  .e.  the  use  of  boron/  epoxy. 

1 1  i  s  know  n  that  the  ampi  i  tude  of  vi  brati  on  i  s  i  nversel  y  proporti  onal  to  the  square  root  of  the 
I  oss  factor.  FI  ence  the  hi  gher  the  I  oss  factor  the  I  ower  the  ampi  i  tude  of  vi  brati  on  or  stress  i  n 
theplate.  Shown  in  Figure3isaplotof  shear  modulusand  I  oss  factor  versus  temperature  for 
a  viscoelastic  damping  material.  Three  distinct  states  exist,  namely,  glassy,  transition  and 
rubbery.  The  maximum  loss  factor  occurs  in  the  transition  state.  It  is  evident  that  these 
vi  scoel  asti  c  dampi  ng  materi  al  s  are  opti  mum  over  a  I  i  mited  temperature  range.  The  response 
i  s  al  so  dependent  on  thefrequency.  FI  ence  the  correct  materi  al  must  be  used  for  the  expected 
environment  (temperature  and  resonant  frequency).  The  material  response  for  different 
damping  materials  is  considered  in  Section  3. 
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Simply  supported  edges 


Clamped  edges 


Figure 2:  Shear  deformation  of  viscoelastic  layer  under  constrained  layer  damping 


Figure  3:  Plot  of  shear  modulus  and  loss  factor  as  a  function  of  temperature  at  constant  frequency  for 
a  viscoelastic  material 
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Tosimplifythe3D  analysis  of  damped  repairs,  an  analytical  solution  to  a  three- layered  plate 
will  be  considered  in  order  to  obtain  the  I oss  factor  and  the  resonant  frequency.  Consider  a 
rectangular  plate  as  shown  in  Figure4with  asingleviscoelasticlayer  and  aconstraining  layer 
with  simply  supported  edges.  Themoduli  of  elasticity  of  the  pi  ate  and  the  constraining  layer 
ar eEi  and  E3  respectively,  and  thedamping  layer  is  characterized  by  its  shear  modulus  G2. 
Hence  the  flexural  stiffness  of  layers  land  3ar  eEih  and  E3I3,  where/i  and  h  are  the  second 
moment  of  area  of  the  cross-sections  of  layer  land  3  respectively. 


View  A  A 

Figure4:  Geometry  for  constrained  layer  damping  of  a  simply  supported  plate 


Let  the  pi  ate  I  ieinthex-y  plane,  then  thegoverning  equation  of  the  displacement,  w ,  inthez 
direction  is  given  by  [p386, 12]: 


d  6w  * 

El  —E  -  g  El , 

dx6 


—  +  m, 

4  1 


d4w 

dx2dt2 


8  w  d2p(x,t )  *  .  . 

=Hnri-* />(*.'> 


where  p(x,t )  is  space-time-dependent  acoustic  pressure,  g*  isthecomplex  shear  parameter 
(defined  later)  and  mt  is  the  total  mass  per  unit  length  of  the  beam.  The  bending  stiffness 
terms  are 

El,  =  EXIX  +  £3/3  (2) 

EIS  =  EIt  (1  +  Y)  (3) 

with  the  geometric  factor  7  [p389, 12]  being 
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Y  _  ji3  E\h\  EJi,  ^ 

Dl(\-v2){Elh]  +E3h3 ) 

The  solution  of  this  equation  leads  to  the  explicit  relationship  for  the  loss  factor  for  a  three 
layer  pi  ate  given  by  [p389, 12], 

n  = _ *e M _  (5) 

,s  i+g„(2  +  i’)  +  ^(i+^)(i  +  r)' 

where  gm  is  the  geometric  shear  parameter  and  rjc  is  the  loss  factor  for  the  viscoelastic 
materi  al .  I  n  ord er  to  accu  ratel  y  p red  i  ct  the  mod al  I  oss  factor,  parti  cu I  arl  y  at  I  ow  temperatu  res 
when  the  loss  factor  of  the  damping  material  is  low,  it  has  been  necessary  to  modify  Mead's 
solution,  [12],  to  include  the  I  oss  factor  of  the  skin  and  the  constraining  layer. 

Denoting  theloss  factors  of  thealuminium,  the  viscoelastic  material,  and  the  boron  patch  as 
rjx,  rj2,  and  r/2 ,  respectively.  TheYoung's  moduli  of  aluminium  (Ex),  the  shear  modulus  of 
the  visco-elastic  layer  {G2 ),  theYoung's  modulus  of  the  boron  patch  (£3 )  are  then  given  by 
£,  =£,(1 +  /?;,)  (6) 

G2  =  G2(l  +  ir/2)  (7) 

E3=E'3(l  +  irj3)  (8) 

with  E\ ,  £’ ,  and  g  denoting  real  parts  of  the  complex  Young's  moduli  of  the  plate,  the 
damping  layer,  and  the  shear  parameter,  respectively.  The  shear  parameter  isthus 
„_G2(\  +  ir/2)(  1  ,  1 


£j  (1  +  //7j  )£j  £  (I  +  ///,)/?, 


=  gO+i?Is) 


where 


T  1  +  ir/n  1  +  ir/2  .  1  +  ir{ j2  1  +  /?;2 

n^=Im  — ; - - —  +  — - - -  /Re  — - - —  +  — - - - 

[_£1(1  +  /'71)/21  £(1  +  ///,)£  J  [£,(!  +  ///, )/?,  £ ( I  +  ///, )/z. 

The  above  expression  leads  to: 

£A(y72  —  ^3  )(1  +  )  +  £3^3  (^2  —  77i  )(1  +  ^3 ) 

£A(!  +  >hJh)(\  +  7i2)  +  +  7i72)(1  +  4) 

It  is  readily  seen  that  the  above  solution  gives  773.  =  ?;2  when  =  ?;3  =  0 . 


The  total  bending  stiffness  of  the  repaired  structured 

El,  =  £,/,  +  £3/3  =  £j(l  +  ///,)/,  +  £3(1  +  />73)/3  =  £/h.(1  +  irj,)  ( 

with 

_  A  A7i  +  £3/3  ?73  ^ 

A/1+A/3 

Themth  complex  frequency,  given  by  Equation  (9.52)  in  ref  [12],  can  bewritten  as 

,v  s"(1+"7s)  ^  ' 

m, 


l  +  g„,(1  +  *75) 
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8m=8  /k}  = 


k]h2 


1  1 

-  +  - 


EA 


ElA  J 


(l-fa 


(15) 


The  modal  loss  factor  rjDC  m  is  the  ratio  between  the  imaginary  part  and  the  real  part  of  the 


expression  on  the  right  hand  side  of  the  equation  (13).Theeffectivewavenumber  kf  will  be 
defined  later.  After  some  simplification,  the  modal  loss  factor  becomes, 


with 


Ui±n 

1  -  7,7  ’ 


(16) 


V  = 


8mVsY 


(17) 


l  +  gm(2  +  7)  +  g^(l  +  7s2)(l  +  7) 

which  is  the  original  solution  by  Mead,  without  considering  the  structural  loss  due  to  the 
aluminium  and  boron  patch. 


Since  77,  «1  for  most  metallic  materials,  the  modal  loss  factor  is  approximately  equal  to 

VnC,m  =  Vt+V  ■  (18) 

Hencethemodal  loss  factor  will  equal  the  loss  factor  duetotheskin  and  constraining  layer.  A 
valueof  77,  for  7075-T6  at  room  temperature  is  avail  able  from  ref.  [30]  where  77  =  25  =0.0106. 
This  will  not  change  significantly  with  temperature  or  frequency.  For  built  up  aircraft 
structures  the  skin  damping  isenhanced  by  frictional  damping  dueto  fasteners  resulting  in 
typi  cal  I oss  factors  of  the  order  of  3%. 

As  indicated  the  geometric  shear  factor  can  be  expressed  as 

g  =  g  /k}  (19) 

and  for  two-dimensional  structures,  the  shear  parameter  is  related  to  the  wave  numbers  in 
two  mutually  perpendicular  directions,  which  is  defined  by 

k2f=k2x+k2y,  (20) 

where  themodal  wave  numbers  inthex  and  ydi  recti  onsforthefundamental  mode  given  by, 
kx  =  ntLx  for  simply  supported,  =  2 n!Lx  for  clamped  (21a) 

ky  =  n!Ly  for  simply  supported,  =  2k IL  for  clamped  (21b) 

with  Lx  and  L  being  the  lengths  of  the  plate  inthex  and  y  directions  respectively. 


It  can  be  shown  that  for  a  three-layer  beam  the  optimum  I  oss  factor  is  given  by  [p392, 12]  and 
is  modified: 


Vs,  max  — 


77 cY 


+  Vt 


(2  +  7)  +  2^/(1  +  7)(1  +  r/2) 

and  an  optimum  valueof  the  shear  parameter  exists,  and  isgiven  by  [p392, 12]: 
1 

8op'  a/(1  +  7)(1  +  T7C2) 


(22) 


(23) 


If  the  design  is  not  close  to  the  optimum  then  consider  a  change  in  7,  in  particular  increase 
the  val  ue  of  fa  and/  or  fa.  There  may  be  an  aerody  nami  c  restri  cti  on  of  the  val  ue  of  fa  such  that 
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the  thickness  of  the  repair  remains  within  the  boundary  layer.  Key  dimensions,  material 
properties  and  parameters  are  shown  in  Table  1.  The  geometric  shear  parameter  gm  in 
equation  [12]  is  a  function  of  the  wave  number  kf ,  which  in  turn  is  a  function  of  the  plate 
dimensions  Lr  and  L  . 

2.1  Design  process 

The  design  process  is  illustrated  in  Figure  5.  Firstly  an  initial  guess  is  required  for  theCLD 
treatment  for  the  structure.  Notethatin  Figure  5  the  sound  pressurelevel  (SPL)  is  required, 
this  will  be  defined  in  Section  4.  The  shear  modulus,  loss  factor,  Young's  modulus, 
temperature,  frequency  and  SPL,  are  required.  From  these  quantities  geometric  shear  factor, 
geometric  factor  and  loss  factor  are  calculated.  Itmay  beuseful  tofirstexaminetheoptimality 
for  this  geometry.  Ifthedesign  isnot  optimum,  thedimensions  of  thedamping  layer  can  be 
adjusted  by  changing  hi,  h2  or  h3.  For  thefirst  design  iteration  neglect  the  optimum.  Firstly 
predi  ct  di  spl  acement  and  frequenci  es.  The  next  step  cal  cul  ates  stress  i  ntensity  dueto  theSPL 
and  residual  thermal  stresses.  Likely  crack  growth  rates  and  reductions  need  to  be  assessed. 
Thedesign  mustbeacceptablethroughouttheentireoperatingtemperaturerange.  Itmay  be 
necessary  to  choose  an  alternativedamping  material  if  results  arenot  acceptable.  Theiterative 
design  process  is  shown  in  Figure  5. 


T ablel:  D  esign  parameters  including  dimensions,  properties  and  parameters 


Physical  dimensions  (Figure 4) 

Material  properties 

Key  parameters 

hi,  h2 ,  h3 ,  hi3 , 

Lx ,  Ly ,  2a 

Ei,E3,G2, 

V1.V3,  A ,  p2  Pi 

Vs  ■  Vc  '  V, ■  gm'  kx'  ky  , 
kf  .  f0  ,  K-rms  ,Y 
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Figure  5: 


Flow  chart  for  design  of  highly  damped  repairs 
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Figure  7:  Dyad  606  material  data,  [14].  (To  obtain  shear  modulus,  multiply  vertical  scale  by  106). 
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Figure8:  Dyad  609  material  data.  (Toobtain  loss  factor  divide  shear  modulus  scale  by  106) 

The  proposed  design  of  the  highly  damped  repair  is  the  configuration  shown 
d ^grammatically  in  Figure6.Thetheory  developed  hereisapplicableto  simply  supported  or 
fully  clamped  rectangular  plates.  However,  for  best  results  the  edge  support  conditions 
should  allow  for  shear  deformation  of  the  viscoelastic  damping  layer,  as  shown  in  Figure  6, 
otherw i  se  the  maxi  mum  possi  bl  e  I  oss  factor  will  not  be  achi  eved .  N  ote  that  both  an  erosi  on 
layer  and  sealant  were  necessary  to  protect  the  repair. 

For  the  design  of  the  repair  the  foil  owing  points  were  considered: 

a)  The  constrai ni ng  layer  had  an  equivalent  stiffness  to  the  plate  being  repaired. 

b)  The  dampi  ng  materi  al  was  chosen  such  that  opti  mal  dampi  ng  occurs  for  operati  on  i  n 
the  expected  temperature  range.  In  this  case  Soundcoat  Dyad  609  was  chosen.  The 
values  of  77,  and  shear  modulus  G2  areobtained  from  Figure7or  Figure8for  Dyad 
606  and  Dyad  609  materials,  respectively  [14],  using  the  procedure  covered  in 
Section  3.  Also  the  thickness  of  the  damping  layer  must  be  specified. 

c)  The  load  to  be  applied  was  the  power  spectral  density  (PSD)  of  the  excitation,  as 
described  in  Section  4. 

d)  Determine  the  computation  of  the  stress  intensity,  stresses,  strains  and  lateral 
deformation  as  described  in  Section  5. 

e)  Assess  the  thermal  residual  stresses  as  presented  in  Section  6. 

f)  Determi  ne  crack  growth  rates  as  detai  I ed  i  n  Secti on  7. 
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2.2  Edge  support  conditions 

Experimental  evidence  indicated  that  fatigue  cracking  usually  occurred  along  the  edges  of  the 
plate  rather  than  in  the  centre  of  the  plate.  This  suggests  that,  as  a  first  approximation, 
clamped  edge  conditions  are  appropriate. 

3.  Damping  materials 

Damping  data  is  shown  in  Figure7and  Figure8for  Soundcoat  Dyad  606 and  609,  [14].  To 
obtain  the  damping  properties  from  these  plots  follow  the  steps,  below  with  reference  to 
Figure  9: 

(a)  For  the  specific  panel  resonant  frequency  (right  hand  scale)  project  horizontally  to 
i  ntersect  the  appropriate  obi  ique  temperature  curve. 

(b)  Now  intersect  the  shear  (storage)  modulus  curve  vertically  and  read  off  the  shear 
Storage  modul  us  from  the  left  hand  scale. 

(c)  Obtain  the  loss  factor  by  intersecting  the  loss  factor  curve  vertically  and  read  off  the 
corresponding  I  oss  factor  from  the  I  eft  hand  scale. 


Ts  T4  T3  T2  Ti 


Frequency 

Hz 


Figure 9:  Procedure  to  determine  damping  properties  at  the  required  frequency  and  temperature. 
H  igher  values  of  thesubscipt  i  for  T.  indicate  increasingly  higher  temperatures 
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4.  Sound  pressure  levels 

Sound  pressure  levels  (SPL)  measured  on  the  external  surface  of  an  F/A-18  inlet  nacelleare 
plotted  in  Figure  10.  The  spectrum  level,  relative  to  the  overall  sound  pressurelevel  (OASPL), 
was  derived  from  in-flight  one-third  octave  SPL  measurements,  [2].  This  spectrum  level  is 
now  used  to  cal  cu  I  ate  the  PSD  of  the  excitation,  OASPL  pressure  levels  of  172. 2d  B 

havebeen  measured  on  theinlet  nacelle[2]. 
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Figure  10:  Spectrum  and  one-third  octave  band  I  e/els  of  sound  pressureover  theexternal  nacelle  inlet 
(whereOASPL=172.2  dB)  [2],  The  relationship  between  the  spectrum  SPL  and  the 
root-mean-square  (r.m.s)  of  the  fluctuating  pressurep,  (Pa),  isgiven  by  [15]  as: 


i  nSPL/ 20-4.65894 
Prms  =  10 


(24) 


and  the  power  spectral  density  of  acoustic  pressure,  i.e.  PSD  of  the  excitation,  at  any  given 
frequency  isgiven  by: 

W0)  =  Pl  =10™"9'3979  (25) 


The  curve  in  Figure  10  has  been  approximated  by  two  straight  lines  defined  by  the  three 
pointslisted  inTable2and  isalso  shown  in  Figure  10.  This  approximate  spectrum  is  used  as 
the  excitation  PSD  for  the  finite  element  analysis. 
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Table2:  Input  power  spectral  density  corresponding  to  an  OA  SPL  of  172.2  dB 


Frequency  (Hz) 

Three  points  to 
approximate  curve 
(dB) 

Pressure  spectrum  level 
(dB) 

Wo) 

(MPa)2/  Hz 

31.5 

-32.2 

172.2-32.2=140 

4.0x10-8 

1000. 

-35.2 

172.2-35.2=137 

2.005x10-8 

8000. 

-48.1 

172.2  -  48.1=124.1 

1.028X1O9 

From  Table  2,  it  is  possible  to  change  the  spectrum  level  SPL  to  any  OASPL  required. 
Validation  studies  outlined  in  Appendix  A  show  that  a  good  match  bdween  experimental 
data  and  predicted  strains  occurs  at  an  OASPL  of  159  dB. 

5.  C omputati on  of  stress  i ntensi ty  factor 


Acoustic  fatigue  is  associated  with  the  random  loading  of  a  structure  over  a  wi  defrequency 
range.  In  the  case  considered  here,  the  loading  is  applied  to  a  rectangular  panel  shown  in 
Figure  4.  The  complete  analytical  solution  of  the  dynamic  response  of  the  panel  is  very 
complex  [19],  However,  inthecaseof  acoustic  fatigue,  anumber  of  assumptions  can  be  made, 
as  the  emphasi  s  here  i  s  to  determi  nethe  crack  growth  dri  vi  ng  force.  Fi  rstly,  theti  me  varyi  ng 
pressurecan  be  considered  to  beuniform  and  in  phaseoverthepanel.  Secondly,  sincethefirst 
mode  contributes  to  98%  of  the  panel  response[3],  the  contribution  of  all  other  modes  will  be 
ignored.  Also  the  variation  of  power  spectral  density  (PSD)  may  be  assumed  to  be  constant 
with  frequency  near  thefirst  mode  frequency.  Furthermore,  itisassumed  that  the  loss  factor 
using  equations  (11-13)  are  applicable  to  a  rectangular  plate  supported  on  all  sides.  The 
solution  [11]  isdirectly  applicableto  simply  supported  edges,  however  a  modification  may  be 
made  to  apply  to  clamped  edges.  Together  with  the  assumptions  of  [20],  the  M  iles  solution 
gives  the  root-mean-square  (in  time  space)  of  the  stress  intensity  response  as: 


K 


rms 


*/&(/„) 


V  2lls 


K 


(26) 


where 


K''ms  is  the  rms  stress  intensity  factor  dueto  the  acoustic  loading, 

K  i  s  th  e  stati  cstressintensityfactorduetotheapplicationofaunitpressure 

on  a  structure, 

S,  (/0 )  is  the  PSD  of  the  excitati  on  at  f0  (Pa 2  /  Hz) , 


/„  is  the  first  resonant  frequency  (Hz)  which  depends  on  boundary 

conditions  and 

rjs  is  the  I oss  factor. 


In  general,  aircraft  skin  panels  are  usually  connected  to  other  plates  or  stiffeners  at  their 
bound ari  es,  consequently  the acousti  c  response  I  i  es  somewhere  between  the  I  i  miti  ng  cases  of 
simply  supported  edges  or  fully-fixed  edges.  Cracking  on  the  edge  of  the  longest  side 
suggests  that  clamped  boundary  conditions  dominate. 


13 


DSTO-RR-0283 


5.1  Un -repaired  plate 

5.1.1  Mode  shape 


1 1  i  s  necessary  to  defi  ne  a  mode  shape  that  sati  sfi  es  the  boundary  conditi  ons  at  the  edge  and 
centre  of  the  plate.  For  clamped  edge  conditions  a  mode  shape  has  been  found  from  [22]  to 
give  good  results.  This  mode  shape  can  be  represented  by: 


w(x,y)  =  w0 


- 

f  \ 

2 

2 

f  ^ 

2 

1-4 

X 

1-4 

UJ 

l4j 

(27) 


w  here  w0  i  s  the  maxi  mu  m  out  of  pi  ane  d  i  sp I  acement  at  the  centre  of  the  p I  ate.  A  ccord  i  ng  to 
M  i  I  es  sol  uti  on,  the  maxi  mum  transverse  di  spl  acement  w0  at  the  centre  of  a  pi  ate  i  s, 


w0 


W(LJLX) 
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EK 


(28) 


where  W(Lyl Lx)  is  given  in  ref  [24];  the  numerical  values  for  a  clamped  panel  can  be 
approximated  by  ref  [35] 


W  =  0.0284e 


-0.754(1,  !LV  y 


Hence  the  displacement  at  any  point  is  given  by, 

w(x,y)  = 


Kw.)  W{yL)  4 


2lts 


4  E\K 


- 

f  \ 

2l 

2 

f  \ 

2 

1-4 

X 

1-4 

y 

v4j 

Computational  validation  is  described  in  TableAl  of  Appendix  A. 


(29) 

(30) 


5.1.2  Resonant  frequency 


For  the  case  of  an  un-repaired  platewith  all  edges  clamped,  itisassumed  that  the  crack  is 
small  enough  that  it  does  not  significantly  affect  the  first  resonant  frequency.  H  encethefirst 
modefrequency  is  given  by: 


in  \ 


D, 


P  h  4 


(31) 


where/?  is  the  density,  Lx  and  Ly  are  the  length  and  width  of  the  plate  respectively.  The 


parameter  7}  =  Eh3  / 12(1-4)  denotesthebendingstiffness.Thegeometryfactor  F(L  /  Lv) 

isgiven  in  ref  [24].  For  the  case  of  clamped  edges,  the  geometry  factor  can  be  approximated  by 
ref  [35], 


F  =  e 


0A71(LXIL  Y- 


(32) 


5.1.3  Strain  and  stress 

The  deformation  of  an  un-repaired  plate  under  acoustic  loading  is  of  pure  bending,  with  the 
bending  strain  along  the  x-di  recti  on  being  given  by, 
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_  \  d  w 
2  dx2 


(33) 


Along  the  clamped  edgey=0,  the  root-mean-square  of  the  bending  strain  becomes 

L\  h  16 


rms  _ 

£  u 


KWo)  r(4}  -V 


2rls 


L.  Exh\  2L\ 


(  T  "j 

2 

12 

A 

-1 

V  Lx  ) 

(34) 


Wi  th  the  strai  n  paral  I  el  to  the  cl  amped  edge  bei  ng  zero,  the  rms  val  ue  of  the  bend  i  ng  stress  i  s 

E, 


_ rms 

u  — 


(l-vf)V  2/7, 


K  \  16 


l:er  2K 


f  x  3 

2 

12 

X 

-1 

l4j 

(35) 


5.1.4  Stress  intensity  factor 


Theanalysisto  determine  the  stress  intensity  factor  for  the  pi  ate  alone  subjected  to  uniform 
pressureis  given  in  ref  [17].  This  analysis  isapplicableto  both  simply  supported  and  clamped 
plates.  The  bendi  ng  stress  i  ntensity  factor  of  a  crack  of  total  length  2 a  is  given  by, 


Kh  =  <jhA7ia 


1  +  v 
3  +  v 


(36) 


Applying  the  M  iles  solution,  the  rms  value  of  the  bending  stress  intensity  factor  becomes, 


jyrms  rn. 

Kb  =  °b 


\na- 


l  +  v 
3  +  v 


(37) 


with  cr™“  being  given  by  equation  (33).  Computational  validation  is  described  in  Appendix 
A. 


5.2  Repaired  plate 


5.2.1  Mode  shape  and  resonant  frequency 


The  same  mode  shape  for  the  un-repai  red  plate  is  used.  For  the  case  of  an  un-repaired  plate 
with  all  edges  clamped,  it  is  assumed  that  the  crack  is  small  enough  that  it  does  not 
si  gni  fi  cantl  y  affect  the  fi  rst  resonant  frequency.  W  hen  the  stiffness  of  the  vi  scoel  asti  c  I  ayer  i  s 
negligible,  i.e.,  zero,  the  first  mode  frequency  is, 


f=hF(L’IL\ 


I^L. 

mA 


(38) 


wherem,  =  plhl  +  p2h2  +  p2h2  isthetotal  mass  per  unit  length,  and  the  bending  stiffness 


D  _  M  i  M 

'  12(  1  -  Vj2)  12(1  -v32) 

and  F(Lyl Ly )  is  given  by  expression  (32)  for  clamped  plate. 


(39) 


The  correct  frequency  for  a  non-zero  or  fi  nite  stiffness  core  can  be  determi  ned  by  the  use  of  a 
non-dimensional  frequency  ratio, 
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Q  =  Jl  +  Y 


g, 


(!  +  gm  (1  +  Us)) 


1  +  2g,n  +  +  7s) 

The  resonant  frequency  for  fi  nite-stiffness  core  i  s  thus 


f°-hF(L’/L\ 


_D, 

m,L l 


-O 


(40) 


(41) 


Results  in  Table  Al,  Appendix  Al,  show  good  agreement  between  solution  (27)  and  the  FE 
analysis. 

5.2.2  Strains  and  stresses  in  a  repaired  plate 


Followingtheabovederivations,  the  root-mean-square  (RMS)  out-of-planedisplacement  w0 
at  the  centre  of  a  repai  red  pi  ate  i  s  gi  ven  by 

1 


wTs  = 


KWoV(^)_ 


L\ 


(42) 


V  2rh  4 '  E$  +  EM:  Q2 

Using  the  same  mode  shape  given  by  equation  (27),  RMS  of  the  out-of-plane  displacement 
w(jc,j>)  can  be  expressed  as, 


wrms(x,y)=wrfl(x)f2(y)  =  w[ 


rms 

M) 


- 

f  \ 

21 

2 

6  A 

2 

1-4 

X 

1-4 

- 

.  L  i 

1  L  , 

_ 

V  x  J 

\  y ) 

(43) 


The  bending  strain  alongy=0can  be  derived  in  a  similar  way  as  for  theun-repaired  plate, 

p\2  rms  j 

nns=_OW  ^ 

dx2  2 

The  r.m.s  membrane  strain  in  the  plate  is  obtained  as  (Mead  and  Markus,  1969) 


rms 

g rms  _  VU  i  _ 


1 


dx 


gdExhx 


mtco2wrms-Dt 


r 

o  w 
dx 4 


+  gD,Y 


^2  rms 

o  w 
dx2 


.(45) 


Since  the  strain  parallel  to  the  clamped  edge  is  zero,  the  RMS  stresses  are 


77  prms 
rms  _  bm 

1-v2 

P  „rms 
rms  _  ■c'l bb 

°b  ~  7  ~2~ 


1  — V, 


(46a) 

(46b) 


Using  strains,  this  expression  is  validated  against  F.E.,  and  the  results  are  given  in  Table  Al, 
Appendix  A. 


5.2.3  Stress  intensity  factor  of  repaired  plate  subjected  to  acoustic  pressure 

The  most  common  typeof  repai  r  appl  i  ed  to  ai  rcraft  ski  ns  areonesi  ded  as  i  1 1  ustrated  i  n  Fi  gure 
11,  si  nee,  i  n  general ,  access  i  s  onl  y  avai  I  abl  eto  the  outsi  de  of  the  ski  n .  A  nal  yti  cal  sol  uti  ons  for 
one-sided  repairs,  known  as theWang-Rose model  areavailablefor  both  in-plane[16, 17]  and 
out-of-plane  loading  [17],  A  significant  feature  of  bonded  repai  rs  is  thatthevariation  of  stress 
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intensity  K(z)  with  crack  length  does  not  increase  indefinitely  with  crack  length  [18,  25]. 
I  nstead,  the  stress  intensity  approaches  an  asymptotic  value.  Thethickness  of  the  adhesive  is 
included  in  the  analysis  to  represent  the  damping  material.  The  computation  of  stress 
intensity  was  based  on  the  work  carried  out  in  ref  [11]. 


Consider  thedamped  plate  shown  in  Figure  11.  Depicted  in  Figure  11a  arethecomponents  of 
thedamped  pi  ate  subjected  to  a  lateral  loading  resulting  in  a  bending  moment.  Intheabsence 
of  a  crack,  the  bending  stress  which  varies  linearly  through  the  pi  ate  thickness  is  shown  in 
Figure  11a.  The  Wang-Rose  crack  bridging  theory  involves  firstly  replacing  a  crack  in  the 
damped  repair  by  springs  along  the  face  and  secondly  use  of  the  lap  joint  analysis  to 
determi  ne  the  I  oad  transfer.  Thi  s  results  i  n  the fol  I owi  ng  equati ons  for  the  membrane  stress 
intensity  Km  and  the  bending  stress  intensity  Kh ,  assuming  the  membrane  and  bending 
stresses  are  constant  over  the  crack  length, 


j^rms  _  /'  rms  r 

v  m  J  mm  m  J  n 


rms 
mb  b 


(47a) 

T^rms  r  rms  .  r  rms  (  AHV\\ 

Kb  =  J  bm  +  Jbb  °b  (47b) 

where  fmm,  fmh ,  fbm  and  fbb  are  the  normalised  stress  intensity  factors.  Plots  of 


fmm  •>  fbm  J  mb  >fbb  are  shown  in  Figure  12-Figure  14  for  temperatures  of  -67  °C,  30  °C  and 
120 °C.  Figure  15  shows  the  variation  of  the  normalised  stress  intensity  factor  with 
temperature.  The  closed  form  sol  uti  ons  arein  close  agreement  with  FE  analysis,  as  shown  in 
TableA2,  Appendix  A. 


In  the  present  case,  as  the  crack  length  increases,  the  prospective  membrane  and  the  bending 
stresses  vary  al  ong  the  crack  I  ength.  I  n  thi  s  case,  a  ful  I  crack  bri  dgi  ng  anal  ysi  s  i  s  requi  red  to 
account  for  the  effect  of  stress  gradient.  For  an  edgecrack  lying  between  (-a ,  L  12  )and  [a  , 


LI  2  ),  the  prospective  membrane  stress  crm  and  the  prospective  bending  stress  ab  both  vary 

with  the  coordinate.  The  crack-bridging  problem,  viz,  the  integral  equations  [17],  can  be 
solved  numerically  using  a  Galerkin  method:  expand  the  unknown  functions  in  terms  of 
Chebyshev  polynomials  and  then  determine  the  coefficients  numerically.  Because  the 
prospective  stresses  vary  with  the  position  al  ong  the  crack  path,  the  matrix  equati  ons  become 
[35] 


i  i 

Ajfj  +  Ajgj  =  —  f  V I  -  f  Uj  (: r)am  (: r)dr  (i,  j  =  0,  l,  2,  ■  ■  ■  N) 

hs  1 

Cjj+D0gJ=±\yli^7uj(r)ab(r)dr  (i,  j  =  0,  1,  2,-N) 

hs  1 


(48) 


where  Uj  denotes  the  Chebyshev  polynomials  of  the  second  kind. 


Figure  16  shows  the  frequency  function  for  the  variation  of  the  first  mode  of  a  rectangular 
plate  for  all  edges  simply  supported  or  all  edges  clamped. 
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'H 


+ 


I  hj  2. 


(a)  2D  view  of  un-cracked  reinforced  case 


Figurell:  Onesided  repair  consisting  of  a  reinforced  beam  subject  to  oneunit  of  lateral  pressure 
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Figurel2:  Resultsfrom  W  ang-Rose  crack  bridging  model  [17]  for  lateral  static  loading  at  -67  °C 


b 

2 


0.24 

0.22 

0.20- 

0.18 

0.16 

0.14 

0.12- 

0.10 

0.08 

0.06- 

0.04- 

0.02 

0.00 


Temp.  =  30°C 
Dyad  609,  4  ply  patch 
■  fmm 
•  fbm 
— A— fmb 
▼  fbb 


■  -  t — ▼ — v- 

7T/ 

V 

- - 


—A — A 


-A - A- 


— I — 

20 


40 


60 


80 


100 


120 


Crack  length  a  (mm) 


Figure  13:  Results  from  the  1/1/ ang-R ose crack  bridging  model  [1 7]  for  lateral  static  loading  at  30  °C 
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Figure  14:  R  esults  from  the  1/1/ ang-R  ose crack  bridging  model  [1 7]  for  lateral  static  loading  at  120  °C 


Figure  15:  Results  from  theW  ang-R  osecrack  bridging  model  [17]  for  lateral  static  loading  and  crack- 
length  a=  25.4  mm 
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Figure  16:  Frequency  function  for  first  mode  of  a  rectangular  plate,  all  edges  clamped  and  all  edges 
simply  supported. 

6.  Residual  thermal  stresses 

Atagiven  operating  temperature,  thethermal  stresses  can  beconsidered  as  the  superposition 
of  thethermal  stresses  due  to  two  steps:  the  curing  of  the  adhesive  (involving  a  localized 
heating-cooling  cycle)  and  a  uniform  temperature  change  from  the  room  temperature  to  the 
operating  temperature.  For  an  ai  rcraft  panel  thatisriveted  on  to  stiff  frames,  the  aircraft  skin 
is  completely  restrained  from  expansion  and  contraction  during  the  heating-cooling  cycle. 
Consequently  no  thermal  stress  will  develop  as  a  result  of  the  bonding  of  the  repair  patch. 

Duri  ng  operati  onstheai  rcraft  ski  n  temperature  can  drop  to  as  I  ow  as  -67  °C.  With  theai  rcraft 
skin  being  stiffened  by  the  repair  and  constrained  by  frames,  a  high  thermal  stress  will  occur 
in  the  repair  and  the  skin.  A  conservative  approach  is  to  assume  that  all  edges  of  the  skin 
panel  are  completed  constrained,  i.e.,  the  total  strain  remains  zero.  The  resulting  thermal 
stresses  in  the  skin  are  given  by, 

°l=°r„=  7 —a,E,(RT  -  T  )  (49) 

1-1/ 

where  RT  and  Toper  denote  respectively  the  room  temperature  and  the  operating 

temperature,  also  thethermal  coefficient  of  expansion  oftheplateisdenoted  by  ax  .Theabove 
expression  furnishes  an  upper-bound  estimate  of  thethermal  stresses. 

A  lower-bound  estimate  of  thethermal  stresses  can  be  obtained  by  assuming  thattheframes 
exp  er i  en  ce  th  e  same  temperatu  re  ch  an  ge  as  th  e  ski  n  s.  I  n  th  i  s  case,  th  e  p  robl  em  can  be  v  i  ew  ed 
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as  a  large  pi  ate  containing  a  smal  I  patch.  Thethermal  stress  in  the  skin,  underneath  the  patch, 

is  given  by  the  Rose  and  Wang  [34], 

_{\  +  vx)(\-aJax)S 

-axhx{Kl  -loper) 


_  T  _ 

C7  —  CT  —  ■ 

xx  yy 


(50) 


2(1-v3)  +  (1  -v\)S 

whereS  denotes  the  stiffness  ratio  given  by, 

S  =  E3hJEJh  (51) 

The  abovethermal  stresses  are  membrane  stresses.  The  stress  i  ntensity  factor  associ  ated  with 
the  thermal  stresses  are 

K  crTxxfmm  (52) 

(53) 

(54) 

The  influence  of  thermal  stresses  on  fatigue  crack  growth  rates  will  be  discussed  in  the 
following  section. 


T 

xx  J  mm 
t^T  _  T  r 
b  xx  J  bm 

Kt  =KTm  +KTb 


7.  C  rack  growth  model 


By  comparing  data  from  standard  crack  growth  laws  [28]  with  acoustic  data  [27]  it  has  been 
found  that  good  agreement  exists.  The  standard  crack  growth  laws  are  intended  to  be  accurate 

from  6  to  50  MPaVm  ,  and  apply  to  in-plane  tension  cycling  only.  The  threshold  stress 
intensity  factor  for  crack  growth  under  high  frequency  out-of-plane  bending  may  be  higher 
[29]  than  for  in-plane  loading.  The  crack  growth  law  for  7075-T6  is  given  by  [29]  as: 


da 

dn 


=  Cf{AK)m  1 


(aO 

2 

l  A K  j 

(mm/  cycle) 


(55) 


where  Cf  =2.1x10  1 ,  m  =  3.02  .  Here  AK  isassumed  to  be  equal  to  A Krms ,  and  AKfh  isthe 
threshold  stress  intensity  factor  for  crack  growth. 

7.1  Crack  growth  threshold 

Data  derived  from  [29]  indicates  that  the  crack  growth  threshold  for  peak  values  of  stress 
intensity,  A Klh ,  factor  is  approximately  10  MPa-Jm  . 
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FigureU:  Effect  of  crack  closure  in  random  bending,  variation  ofK  with  time  [30] 

7.2  Crack  closure 

It  is  known  that  under  cyclic  lateral  loading,  plasticity  will  occur  at  the  crack  tip  [30].  Asa 
result,  during  thesecond  tension  cycle  there  will  bea  del  ay  until  the  crack  opens.  Thefracti  on 
oftheload  atwhich  the  crack  opens  is  defined  as  (3  .  Asshown  in  Figure  17,  the  time  hi  story 
for  a  random  event  in  which  a  maximum  peak  occurs  is  defined  by  a  Krms .  From  [30]  the 
effective  (peak)  stress  intensity  factor  dueto  crack  closure  is  given  by: 

AKeff  =  Kp  -  a  (3  Krms  (56) 

Using  the  values  of  a  and  [3  availablefrom  [30],  i.e.,  a  =3.75and  [3  =0.3  respectively,  the 
effective K  is 

AKeff=Kp  -  1.125  Krms  (57) 

7.3  Crack  growth  for  narrow  random  crack  growth 

1 1  i  s  al  so  kn  ow  n  that  d  u  r  i  n  g  cy  cl  i  c  I  ateral  ben  ding  crack  cl  osu  r  e  of  th  e  p  I  ate  may  occu  r  w  h  en 
thesurfaceof  the  plate  is  in  tension.  This  results  from  theplasticzoneatthetip  of  thegrowing 
crack.  Since  thisplasticzonewill  remain  in  the  wake  of  the  crack,  unloading  of  the  crack  will 
result  in  crack  closure  before  complete  unloading  occurs.  From  [30]  a  peak  probability 
distribution  has  been  used  as  representative  of  acoustic  loading  and  exhibits  a  "truncated" 
probability  distribution  in  which  peaks  above  a  certai  n  level  do  not  occur.  From  flight  trials 
data  [31]  the  ratio  of  peak  val  ues  to  rms  val  ues  is  5. 
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7.3.1  Case  A 


If  the  stress  i  intensity  factor  at  which  the  crack  opens  ( aj3Krms )  is  greater  than  the  peak  val  ue 

of  thethreshold  stress  intensity  V2  K'rms  then  peaks  I  ess  than  a[iKrms  will  not  contribute  to  the 
crack  growth. 


i.e  aPK,m>4lK, 


(58) 


da  _  C  a2p2  a 2  a2 J32 

~dN~  (\.-p2y  rms  2.  ’T’  2. 


(59) 


7.3.2  CaseB 


If  a/3Knns  is  less  than  V2  K‘rms  then  this  will  not  result  in  crack  growth,  dueto  a  threshold 
effect.  In  this  case, 


da 

~dN 

where 


C 


2  o2 


(i-AT 


-k: 


r  ,  or  B  ^ 

1  - 


a 

2. 


(K,,)2 

K: . 


(60) 


h 

I  {n,  d,  h,  l  )  =  | (  tU2  -  d112)"  e ^  dt  . 

i 


(61) 


The  alternative  to  solving  these  equations  is  to  design  the  repair  such  that  K  <  10  Kp ,  for 
7075-T6. 


8.  Acoustic  loading  on  F/A-18 

Whiledataexist[2]fortheSPL  on  the  inlet  nacelle  region,  there  is  no  data  for  the  aft  fusel  age 
region.  It  was  reported  that  the  airframe[2]  was  designed  for  an  overall  sound  pressurelevel 
(O  A  SP  L )  of  160  dB.FIighttrials[31]inwhichadosi  meter  w  as  i  n  stal  I  ed  to  measu  re  th  e  strai  n  s 
on  two  adjacent  skin  panelswith  and  without  repair,  respectively,  have  provided  both  peak 
and  r.m.s.  strain  data  on  the  aft  fusel  age  region.  For  the  un- repaired  panel  two  strain  gauges 
were  used,  with  one  gauge  at  the  centre  of  the  panel  and  the  other  at  the  edge  of  the  panel . 
Sincethecentregaugewasin  aregion  of  low  strain  gradient,  the  measured  strai  ns  are  more 
reliable  than  the  gauge  close  to  the  edge.  Assuming  that  the  mode  shape  is  the  fundamental 
mode,  it  is  possibleto  back-calculate  the  OASPL  from  the  centre  strain  gauge  data.  These 
calculations  aredescribed  in  Appendix  A  and  show  thattheOASPL  in  the  region  of  the  aft 
fusel  age  is  159  dB,  very  close  to  the  design  level  of  160  dB. 
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9.  Highly  Damped  Repair,  Example 

A  n  examp  I  e  probl  em  su  ch  as  a  h  i  gh  I  y  d  amped  repai  r  for  the  aft  f  usel  age  of  the  F/  A  - 18  w  i  1 1  be 
consi  d  ered .  D  ata  need  ed  for  th  i  s  d  esi  gn  are  show  n  i  n  Tabl  e  3,  Tabl  e  4  and  Tabl  e  5.  The  vari  ou  s 
parameters  used  i  n  the  anal  ysi  s  are  show  n  i  n  Tabl  e  7.  The  expected  operati  ng  temperatu  re  of 
theF/A-18inthevicinityoftheaftfuselageis-67to+1200C. 


T able  3:  D  etails  of  parameters  for  the  highly  damped  repair  example 


Item 

Value 

Description 

AT 

40  °C  * 

From  curing  temperature  to  operating  temperature. 

ax 

23.X106/  C 

Thermal  coefficient  of  expansion  of  plate. 

a3 

5.57x106/  =c 

Thermal  coefficient  of  expansion  of  constraining  layer. 

K 

0.9  mm 

Thickness  of  plate. 

h2 

0.508  mm 

Thickness  of  clamping  layer. 

K 

0.508  mm 

Thickness  of  constraining  layer. 

Ex 

71000  M  Pa 

Young's  modulus  plate. 

E, 

113300  M  Pa 

Modulus  of  constraining  layer. 

0.3 

Poisson's  ratio  for  plate. 

^3 

0.0535 

Poisson's  ratio  for  constraining  layer. 

Pi 

2.77x109  kkg/  mm3 

Density  of  plate. 

Pi 

2.03x109  |<kg/  mm3 

Density  of  constraining  layer. 

Pi 

1.10x109  kkg/  mm3 

Density  of  viscoelastic  material. 

L  x 

175  mm 

Length  of  plate. 

Tv 

176  mm 

Width  of  plate. 

\K0/cr\ 

Coefficients  from  Figures  12-15,  depending  on  crack  size. 

m 

3.02 

Parameter  for  crack  growth  law  for  7075-T6. 

cf 

1.37x105 

Parameter  for  crack  growth  law  for  7075-T6. 

Kf 

63.9 

Parameter  for  crack  growth  law  for  7075-T6. 

di 

1.0* 

Loss  factor  for  viscoelastic  damping  layer. 

Vi 

0.013 

Loss  factor  for  plate. 

% 

0.004 

Loss  factor  for  composite. 

g2 

10  MPa* 

Shear  modulusfor  damping  layer. 

*  vari  abl  es  for  exampl  e  probl  em 
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Table  4:  L  oss  factor  and  shear  modulus  for  Dyad  606  versus  operating  temperature  corresponding  to 
cure  temperature  of  90  °C,  from  Figure  7 


Operating  temp.  °C 

Temperature  Change  °C 

Loss  factor  77  +  77, 

G2  (MPa) 

-67 

157 

0.01912 

300 

-40 

130 

0.01912 

300 

-20 

110 

0.01912 

300 

0 

90 

0.01912 

300 

20 

70 

0.02025 

300 

30 

60 

0.1749 

70 

40 

50 

0.4171 

17 

50 

40 

0.3619 

5 

60 

30 

0.1385 

2 

70 

20 

0.04151 

0.7 

80. 

10 

0.02391 

0.4 

90 

0 

0.02146 

0.4 

Table  5:  L  oss  factor  and  shear  modulus  for  Dyad  609  versus  operating  temperature  corresponding  to 
cure  temperature  of  90  °C,  from  Figure  8 


Operating  temp.  °C 

Temperature  Change  °C 

Loss  factor  r/  +  rjt 

G2  (MPa) 

-40 

130 

0.02109 

300 

-20 

110 

0.02109 

300 

0 

90 

0.02109 

300 

20 

70 

0.03224 

300 

30 

60 

0.04088 

280 

40 

50 

0.05717 

200 

50 

40 

0.09915 

130 

60 

30 

0.14563 

85 

70 

20 

0.23625 

50 

80. 

10 

0.36454 

25 

90 

0 

0.40449 

5 

T able  6:  PSD  of  the  excitation  for  plate  alone,  also  the  repaired  case,  interpolated  from  Figure  10 
from  results  for  inlet  naceile[2] 


OASPL  (dB) 

S,  (197  Hz) 

Sj  (344  Hz) 

159 

2.196x10-9 

2.015x10-9 

26 


DSTO-RR-0283 


T able  7:  T erminol ogi es  used  in  M  athematica  program 


Quantity 

Equation/comment 

s 

47,  Stiffness  ratio 

Q 

37,  Frequency  correction  for  non-zero  core 

fo 

29,  Natural  frequency  for  plate 

fo 

35  N  atural  frequency  for  repai  red  plate 

PSD 

24,  Power  spectral  density 

8m 

18,  Geometric  shear  factor 

k) 

19,  Wavenumber 

Y 

3,  Geometric  factor 

11, 

12,  Loss  factor  due  to  skin  and  patch 

11c 

Loss  factor  for  viscoelastic  material  see  Figures  7  or  8. 

lls 

10,  Loss  factor  due  to  ,  r/2 ,  //, 

Vdcw 

15,  overall  loss  factor 

§  opt 

22,  Can  be  ignored,  this  is  only  needed  for  optimum  design 

Vs  max 

21,  Can  be  ignored,  this  is  only  needed  for  optimum  design 

% 

28,  42  transverse  displacement 

rms 
^  m 

43,  Membrane  strain  in  patched  plate 

rms 

bb 

42,  Bending  strain  in  patched  plate 

_ rms 

am 

45a,  M  embrane stress  in  patched  plate 

rms 

45b,  Bending  stress  in  patched  plate 

l s  rms 

Km 

46a  Membrane  stress  intensity  factor  for  repaired  panel 

jy-  rms 

Kb 

46b  Bending  stress  intensity  factor  for  repaired  panel 

^ , 

69,  Delta,  peak  stress  intensity  factor  (ratio  peak/  rms  is  5) 

The  analyses  described  in  Section  5  have  been  implemented  using  a  symbolic  solver 
Mathematica.  Consider  a  plate  of  dimensions  Lx  =  176,  L  =  175  mm  and  thickness 

hx  =  0.9  mm  with  a  repair  of  Lx  =  176,  L  =  175  mm  and  thickness  h3  =  0.508  mm,  which 

corresponds  to  fou  r  I  ayers  of  boron  for  the  constrai  ni  ng  I  ayer  as  show  n  i  n  Fi  gure4.  The  ai  m  i  s 
to  achieve  a  neutral  axis  within  the  viscoelastic  material  in  order  to  maximise  the  shear, 
although  distances  away  from  the  neutral  axis  are  not  far  from  the  optimum.  Consider  an 
operating  temperature  of  -40  to  90  °C.  In  this  case,  Dyad  609  is  used  due  to  durability 
concerns  of  Dyad  606  at  high  temperature.  Dyad  609  material  only  comes  in  at  thicknesses  of 
0.508  and  1.106  mm.  A  thickness  of  0.508  will  be  used.  For  the  boron  doubler,  a  stacking 
sequence  of  [0/  90/  0/  90]  will  be  used,  resulting  in  a  laminate  modulus  of  114  GPa 

Firstly,  from  experimental  data  [29],  both  the  loss  factor  and  the  resonant  frequency  change 
with  increasing  moisture  content.  N  ote  that  the  un-conditioned  specimen  gives  higher  loss 
factors  than  data  provided  by  the  manufacturer.  The  difference  is  due  to  moisture  uptake 
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whereas  the  manufacturer's  data  is  for  low  humidity.  A  change  of  resonant  frequency 
indicated  that  the  shear  modulus  also  changed  with  moisture.  These  are  factors  that  may  have 
to  be  considered  in  determining  the  real  loss  factor.  For  this  exercise  the  manufacturers  data 
will  be  used,  see  Figure  18. 


o 

TO 

4— 

CO 

CO 

o 


20  30  40  50  60  70  80 


Operating  temperature  °C 


Figure  18:  L  oss  factor  versus  operating  temperature  atOASPL  =159  dB,  for  manufacturers  data,  and 
conditioned  data,  this  corresponds  to  a  4  ply  configuration 

The  equations,  with  parameters  given  in  Table  7,  can  be  solved  numerically  or  symbolically 
using  Mathematica™ .  The  stress  intensity  factors  can  then  be  computed  for  a  range  of 
operating  temperatures  and  SPL's.  The  corresponding  PSD  of  the  excitation  istabulated  in 
Table6foraplatewith  and  without  repair.  The  material  properties  of  Dyad  606  and  609  are 
listed  in  Table  4  and  Table  5  respectively,  for  various  operating  temperatures.  These  results 
are  graphically  displayed  in  Figure  19  and  Figure  20  for  OASPL  of  159  dB  and  operating 
temperatures  from  -67  to  120  °C.  In  Figure  20  the  peak  microstrain  for  Dyad  606  and  609  is 
2000  and  955  respectively.  Note  that  the  strain  level  of  1750  microstrain  corresponds  to  the  un¬ 
patched  panel  derived  from  flight  trials  results.  Cleary  the  Dyad  609  patch  has  reduced  the 
strai  n  by  a  factor  of  al  most  2.  N  ote  that  the  Dyad  606  desi  gn  i  ncreases  the  maxi  mum  strai  n  at 
high  temperatu res  to  a  I  ev el  exceed i  ng  thefl  i ght  tri al  strai  ns  on  the  u n- patched  panel .  Thi s  i  s 
the  result  of  the  damping  material  having  a  very  low  shear  modulus  at  higher  temperatures, 
with  the  additional  mass  of  the  damped  layer  increasing  resonant  response.  Also  shown  in 
Fi gure  19  i s  thethreshol d  val  uefor  crack  growth.  A s  di scussed  Secti on  7.1,  thethreshol d  val  ue 
of  10  MPa  m0-5  was  determined  from  an  experimental  study,  ref  [29].  It  is  evident  that  from 
Fi  gure  21  that  the  4  pi  y  Dyad  606  desi  gn  w  i  1 1  not  survi  ve  the  hi  gh  temperatu  re  I  oad  i  ng.  A I  so 
the  Dyad  609,  4  ply  design,  is  just  above  the  threshold  at  low  temperatures.  Flowever,  it  is 
expected  thatwhen  using  Dyad  609asix-ply  design  will  suffice,  see  Section  9.6 
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Figure  19:  Predicted  rms  stress  intensity  factor  versus  operating  temperature  at  OASPL=159  dB, 
using  Dyad  609  or  606,  for  crack  length  2a=50  mm,  this  corresponds  to  a  4  piy 
configuration  wherethermal  stresses dueto uniform  temperaturechangehavebeen  ignored. 

9.1  Adhesives  for  the  repair 

Further  considerations  concerning  thereal  operating  temperatures  were  made  in  [33].  For  the 
col dest  day  the  enti  re  ai  rcraft  ski  n  temperature  i s  -67  °C .  For  the  hottest  day  of  the  year  the 
repair  was  adjacent  to  both  103  °C  and  141°C  temperature  regions.  As  a  result  an  average 
hottest  day  was  taken  to  be  122  °C,  and  the  analysis  was  carried  out  at  120  °C.  The  adhesive 
used  to  construct  the  patch  was  FM  300,  which  was  co-cured  with  the  boron  and  damping 
material  at  177  °C.  The  adhesive  used  to  bond  the  repair  to  the  aircraft  skin  was  FM  300-2  si  nee 
this  adhesive  would  perform  within  the  range  of -67  °C  to  120  °C.  Figure  24  and  Figure  25 
show  schematically  thestructural  components  of  the  patch  system,  and  thethicknessof  each 
component.  In  this  section  adhesive  shear  stresses  due  to  an  acoustic  loading  will  be 
evaluated. 
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Plies:  4,  2a=50mm, 


-80  -60  -40  -20  0  20  40  60  80  100  120  140 
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Figure20:  Predicted  peak  microstrain  for  repaired  panel  versus  operating  temperature  for 
OASPL=159dB,  using  either  Dyad  609  or  Dyad  606.  This  corresponds  to  4  ply 
configuration. 


Figure21:  Through  thickness  view  oftheindividual  structural  components  for  therepair  oftheplate 
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The  dimensions  and  material  properties  are  shown  in  Figure8and  Figured 


Boron 

FM300 

h3 

,  ^  a2 

Dyad  609 

h2 

FM300 

FM300-2 

Skin 

^a2 

Figure22:  D  iagram  showing  thicknesses  of  patch  components 


T abl  e  8:  P  arameters  for  adhesi  ve  shear  stress 


Temp. 

°C 

Gul 

(MPa) 

FM  300-2 

Gfl2 

(MPa) 

FM  300 

g2 

(MPa) 

tal 

(mm) 

t  9 

a  2 

(mm) 

h2 

(mm) 

120 

199 

199 

1 

0.127 

0.127 

0.508 

0 

415 

4150 

300 

0.127 

0.127 

0.508 

Table9:  Adhesives  shear  stress  dueto  acoustic  loading,  using  parameters  from  Table8 


Ex 

(MPa) 

E, 

(MPa) 

K 

(mm) 

(mm) 

°  m 

(MPa) 

°b 

(MPa) 

T  * 
ip 

(MPa) 

**T 

P  peak 

(MPa) 

71000 

114000 

0.9 

0.762 

0.544 

14.662 

14. 

0.304 

71000 

114000 

0.9 

0.762 

6.621 

6.036 

43. 

3.540 

*  i  s  the  pi  asti  c  val  ue  of  the  shear  stress  for  the  ad hesi  ve 

**  is  the  peak  value,  rpeak  =5  *t(  Mathematica  :D  amped  Repair  16e_609b.nb) 


The  shear  stress  in  the  adhesive  due  to  acoustic  loading,  which  induces  membrane  and 
bending  stresses  ( cr„,  and  crb ),  is  given  by, 


Gil 


T  =  ■ 


cr,„  +  cr. 


where 


K 

g: 


tal  +2-^- 


\Ga  i 


’a  2 


+ 


'2  J 


(62) 


(63) 


/  = 


tg\  ^al 


h2  ^  ExE2hxh2 
G2){E  \hx  +E3h3) 


(64) 
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Theparameters  GaX,Gal  and  G2  denote  theadhesiveshear  moduli  for  FM  300-2,  FM300,  and 
the  damping  layer  respectively.  The  Young's  moduli  of  the  skin  and  the  patch  are  denoted 
respectively  as  i^and  E2.  Parameters  hx,  h2  and  h2  denote  the  thickness  of  the  skin,  the 
damping  layer,  and  the  patch  respectively.  The  thickness  of  the  FM 300-2  and  FM300 
adhesives  are  taX  and  ta2.  Parameters  aAlum  and  aBoron  denote  the  coefficients  of  thermal 
expansion  for  the  skin  and  patch  respectively. 

For  acoustic  loading  the  shear  modul  us  of  theadhesive  varies  with  temperature,  seeFigure23 
for  the  stress  strai  n  curve. 


Figure23:  Stress  strain  diagram  for  FM  300  for  the  temperature  range  shown 


Cal  cul  ati  ons  i  n  Tabl  e  8  and  Tabl  e  9for  acousti  c  shear  stresses  i  n  the  ad  hesi  ve  have  show  n  that 
at  an  operati  ng  temperature  of  0  or  120  °C,  both  FM  300  and  FM  300-2  areoperati  ng  withi  n  the 
elastic  region,  hence  no  fatigue  damage  will  occur,  see  Table  8  and  Table  9.  The  FM300  is 
cured  in  the  autoclave  at  177  °C,  together  with  the  boron  and  damping  material.  FM300-2is 
used  to  bond  the  patch  to  the  ai  rcraft  ski  n  at  a  temperature  up  to  120  °C . 
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9.2  Alternating  stress  intensity  factor 

The  effect  of  low  temperatures  may  result  in  tensile  residual  stresses  which  would  keep  the 
crack  open,  thus  increasing  crack  growth  rates.  In  this  case  a  simple  and  conservative 
approach  i  s  to  assume  that  negative  stress  i  ntensity  factor  does  not  contri  buteto  crack  growth. 
In  this  case,  the  effective  crack  growth  driving  force  can  be  expressed  as 

A Keff  =  min(Xmax ,  Kmax  -  Kmm  )  ,  Kmm  >0  (65) 

Under  acoustic  loading,  theminimumstressintensityfactor  Kmin  isequal  in  magnitudetothe 
maximum  stress  intensity  factor  but  opposite  in  sign,  i.e.,  Kmin  =  -Kmax  .  In  the  presence  of 
thermal  residual  stress  intensity  factor  KT ,  the  effective  range  of  stress  intensity  factor  is: 

AKeff=2Kmax,  KT>Kmax  (66a) 

AK<r=Kma+KT,  KT<Kmax  (66b) 

A  conservative  assumption  for  the  KT  <  Kmax  case  i s: 

A Keff  =  Kmax  (67) 

9.2.1  Small  patch 

When  asmall  patch  isapplied,  duetothelack  of  constraint,  high  residual  thermal  stresswill 
exi  st  in  the  skin  after  the  bonding  of  a  Boron  patch.  As  shown  in  Table  10,  the  thermal  stress 
intensity  factor  is  approximately  equal  to  the  peak  stress  intensity  factor  at  20  °C .  Therefore, 
the  effecti  ve  range  of  the  stress  i  ntensity  i  s: 

A Keff  =  Kmax  +Kr*  2 Kmax  (68) 

This  means  that  small  patcheswill  have  large  thermal  residual  stresses  and  may  not  prevent 
crack  growth. 

9.2.2  Large  patch 

For  a  large  patch  that  covers  the  enti  re  bay  (panel )  that  i  s  riveted  to  the  ai  rframe,  no  thermal 
expansion  or  contraction  would  occur  during  the  bonding  process.  Asa  result,  the  repaired 
skin  will  not  experience  any  thermal  residual  stress,  i.e,  KT  =  0  .  For  proof  seeAppendixA.3, 
in  this  case: 

A Kp  =  Kmax  (69) 

Itisdear  that  the  effective  crack  growth  driving  forcewould  be  much  less  than  that  pertinent 
to  a  smal  I  patch.  Thi  s  i  s  consistent  with  the  previ  ous  fai  I  ed-attempt  of  repai  ri  ng  the  acoustic 
fatigue  crack  with  asmall  Boron  patch  that  covered  only  the  cracked  region. 
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9.2.3  Case  KT  >  Kmm 

At  low  temperatures  thermal  mismatch  may  result  in  a  high  stress  intensity  factor  for  the 
combination  of  aluminium  and  boron.  As  result  during  each  cycle,  AKeff=2Kw. 

9.2.4  Case  KT  <  Kmia 

Inthiscaseat  high  temperatures  thevalueof  KT  may  be  negative.  Sincecrack  closurecannot 
occur  then,  Kr  =  0  and  A Keff  =  Kmix . 

9.3  Repair  for  case  of  uniform  change  of  temperature 

Consider  the  case  for  an  aircraft  which  is  not  subjected  to  acoustic  loading  or  flight  loads  and 
does  not  contain  a  patch.  Now  consider  two  points  on  the  structure  located  far  apart  at 
oppositeendsofthestructure.  Uniform  heating  or  cooling  of  thewholestructurewi  1 1  result  in 
a  rel  ati  ve  i  ncrement  i  n  d i  spl  acement  between  thetwo  poi  nts.  I  f  these  poi  nts  are  brought  near 
each  other  then  the  relative  displacement  will  decease.  At  the  limit  at  which  no  separation 
exists,  the  relative  displacement  is  zero.  If  a  crack  is  introduced  into  the  structure  and  each 
point  is  located  on  the  opposite  crack  face  then  heating  or  cooling  will  not  result  in  any 
rel  ati  ve  d  i  sp  I  acement  of  the  crack  faces. 

Ifwenowapplya  patch  whi  ch  i  s  smal  I  i  n  compari  son  to  the  comp  I  ete  ai  rcraft  structure  then 
the  sti  ff  n  ess  of  th  e  str  u  ctu  re  w  i  1 1  not  chan  ge.  A  s  a  resu  1 1  th  e  crack  openingdispl  acement  will 
be  almost  zero.  Hence  uniform  heating  or  cooling  will  not  affect  the  crack.  The  only  crack 
opening  will  occur  from  acoustic  excitation.  Flight  loads  are  small  in  comparison  to  high 
acoustic  loading. 

9.4  Residual  stresses  from  bonding  the  patch 

Bonded  repai  rs  i  nvol  veheati  ng  thepatch  and  ski  n  up  to  temperatures  of  80-120  °C  for  8  hours 
until  the  patch  is  cured.  Due  to  thermal  mismatch  between  the  boron  and  skin  there  is 
potential  for  residual  stress.  Consider  the  specific  case  of  a  patch  applied  to  a  complete  panel 
which  issurrounded  by  structure  consisted  of  frames  and  sheet  stringers  bolted  together.  In 
thisparticular  case  i  t  can  be  sh  ow  n  ( see  secti  on  A .  5  i  n  the  A  p  pend  i  x  A )  th  at  th  e  resi  d  u  al  stress 
is  zero. 

9.5  0  perati  ng  tern peratu  res 

As  menti  oned  i  n  secti  on  9. 1  the  maxi  mum  and  mi  ni  mum  temperatures  i  n  thi  s  area  are -67  °C 
and  120  °C  respectively.  The  -67  °C  temperature  corresponds  to  the  lowest  airspeed  and 
engine  power  that  will  maintain  the  aircraft  in  straight  and  level  flight  at  Flight  Level  (FL)500. 

I  n  this  condition  the  manoeuvre  capabi  I  ity  is  limited.  The  maxi  mum  temper  atu  re  corresponds 
to  Mach  2atFL  350.  Whilethemanoeuvrecapability  for  the  coldest  temperatureis  limited,  it 
i  s  n  ot  n  ecessary  to  u  se  a  I  oad  al  I  ev  i  ati  on  factor .  N  ote  th  at  th  e  rel  ati  on  sh  i  p  of  sh  ear  mod  u  I  u  s  to 
temperatu  re  for  F M  300  i  s: 
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G  =  415-1. 8T  (70) 

Thedamping  material,  Dyad  609,  ismost  effectiveover  the  range  of  temperature  20- 120  °C.  At 
I  ow  temperatu  res  the  I  oss  factor  of  the  repai  r  w  i  1 1  on  I  y  be  0. 004.  A  vai  I  abl  e  d  ata  i  nd  i  cates  that 
the  properties  of  Dyad  609  do  not  change  at  low  temperatures.  Hence  at  low  temperatures  the 
repair  is  not  damped.  The  patch  works  as  a  'stand  off'  repair  at  low  temperatures.  The 
thi  ckness  of  the  d  ampi  ng  I  ayer  i  ncreases  the  bend  i  ng  stiffness  suffi  ci  entl  y  to  reduce  the  stress 
intensity. 

9.6  Results 

Firstly  consider  an  un-damped  patch  in  Figure24.  This  is  a  classic  examplethat  occurs  when 
un-damped  patches  are  used.  Onerepair  is  a  simple  un-damped  boron  patch  repairwhilethe 
other  is  an  un-damped  aluminium  patch.  Both  repairs  result  in  stress  intensity  factors  that  will 
result  in  crack  growth.  The  aluminium  pi  ate  is  bonded  to  the  skin,  but  the  adhesive  softens 
with  rise  of  operating  temperature.  Almost  the  same  result  will  occur  with  a  bolted  plate 
although  the  softening  will  not  occur,  however  the  end  result  is  much  the  same. 
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Figure24:  Predicted  peak  stress  intensity  factor  versus  operating  temperatu  re  for  un-damped  boron 
repair  ,8  piies,  and  un-damped  aluminium  plate  repair 

The  results  for  the  acoustic  loading  are  shown  in  Figure  25  and  indicate  the  performance  of 
the  4,  6  and  8  layer  designs.  Note  that  the  F.E.  results  for  the  4  layer  patch  show  good 
agreement  with  the  analytic  values.  Clearly  both  the  6  and  8  ply  designs  survivethis  loading. 
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Acoustic  loading:  OASPL  =  159  dB 


Temperature  °C 


Figure25:  Predicted  peak  stress  intensity  factor  versus  operating  temperature  at  OA  SPL  =159  dB, 
using  Dyad  609,  for  crack  length  2a=50  mm.  This  corresponds  to  4,  6  and  8  ply 
configuration. 


Acoustic  loading:  OASPL  =  159dB 


Temperature  °C 


Figure26:  Predicted  peak  stress  intensity  factor  versus  operating  temperature  at  OA  SPL  =159  dB, 
using  Dyad  609,  for  crack  length  2a=75  mm.  This  corresponds  to  4,  6  and  8  ply 
configuration. 

It  is  known  that  as  the  crack  grows  the  stress  intensity  in  a  panel,  surrounded  by  structure, 
will  decrease  as  the  crack  tip  approaches  the  boundary,  see  Figure  A6  in  Appendix  A.  The 
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analysis  is  capable  of  predicting  the  change  in  stress  intensity  as  the  crack  grows.  This  set  up 
in  Correction  Tables,  and  also  requires  a  Fortran  program  to  compute  the  new  stress  intensity. 

For  a  crack  sizeof  2a=75  mm,  the  results  areshown  in  Figure  26.  All  threedesignsfor4, 6  and 
8  layers  survive  this  loading.  As  the  crack  length  extends  to  2a=100  mm  in  Figure  27,  and 
2a=150mm  in  Figure  28  all  designs  survive.  It  is  evident  that  the  6  layer  is  always  below  the 
threshol  d  for  crack  growth.  FI  ence  a  6 1  ayer  desi  gn  i  s  practi  cal  si  nee  it  i  s  capabl  e  of  stoppi  ng 
very  large  cracks,  with  no  crack  growth,  and  hence  has  an  infinite  fatigue  life. 


-100  -50  0  50  100  150 

Temperature  °C 


Figure27:  Predicted  peak  stress  intensity  factor  versus  operating  temperature  at  OA  SPL  =159  dB, 
using  Dyad  609,  for  crack  length  2a=100  mm.  This  corresponds  to  4,  6  and  8  ply 
configuration. 


Temperature  °C 


Figure28:  P  redicted  peak  stress  i  ntensity  factor  versus  operating  temperature  for  damped  repairs  for  4, 
6  and  8  ply  lay-up.  The  crack  length  is2a=150  mm. 
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O  veral  I  the  resu  I  ts  show  that  the  anal  yti  cal  sol  uti  on  u  nd  eresti  mates  the  F  .E .  resu  I  ts  si  i  ghtl  y  at 
low  temperatures  and  a  safety  factor  1.05  is  needed.  At  high  temperatures,  a  factor  of  1.2  is 
required. 

Data  for  all  plotted  results  on  Figures  25- 28  arecontained  in  AppendixA,TablesA10-24.  All 
equations  in  this  report  have  been  solved  using  Mathematica. 


Table  10:  Configurations  of  repair 


Mathematica  file 

Number  of  layers 

Lay  up 

Damped  Repair  609a.nb 

4 

0,90,90,0 

Damped  Repair  609a.nb 

6 

0,90,0,90,0,90 

Damped  Repair  609a.nb 

8 

0,90,0,90,90,0,90,0 

Shown  in  Figure  29  is  the  predicted  variation  in  rms  strain  in  the  plate  versus  operating 
temperatures  for  a  4,  6  and  8  ply  configuration. 


OASPL =  159  dB 


Figure 29:  Residual  predicted  rms  strain  versus  operating  temperature  at  OASPL=159  dB,  using 
Dyad  609,  for  crack  length  2a=50  mm.  This  corresponds  to  4,  6  and  8  piy  configuration. 

9.7  Strength  of  Dyad  609 

Manufacturer's  data  for  Sound  Coat  Dyad  609  gives  the  tensile  strength  as  17.2MPa.  Dueto 
the  envi  ronmental  affects  a  knock  down  factor  of  60  i  s  used .  The  repai  r  i  s  I  ocated  on  the  si  de 
of  the  fusel  age  and  the  edge  is  protected  from  erosion  by  sealant.  Also  the  repair  is  located 
within  theboundary  layer.  As  a  result,  the  largest  suction  that  can  be  achieved  is  in  a  vacuum, 
i.e.  0.097  MPa.  Flence  the  margin  of  safety  is  17.2/(60x0.097)-!  =  1.96 
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The  only  other  source  of  loading  would  come  from  a  sideslip  resulting  in  compression 
loading.  Assume  a  side-slipeangleof  45  degrees,  the  load  for  a  speed  of  700  m/  s  and  for  a 
patch  of  size  303x305  mm  is  0.5  p  V2  =  0.2\2MPa  .The  margin  of  safety  against  compression 

is  17.2/(60x0.213)  -  1  =  0.35. 

9.8  Cost  savings  using  highly  damped  repairs 

In  Appendix  B  costs  are  given  for  the  standard  RAAF  scab  metal  repair  and  the  highly 
damped  repai  r.  These  costs  show  that  w  hi  I  e  damped  repai  r  was  si  i  ghtl  y  I  ower  i  n  costs,  I  ay-up 
oftheboron/  epoxy/  Dyad609  repair  contributes  to  thelabour  costs.  Mass  production  would 
result  in  cost  savings  allowing  patches  to  be  cut  out  from  large  sheets  of  patch  material.  The 
main  difference  between  the  repai  rs  is  fatigue  life.  It  has  been  shown  in  this  report  that  with 
appropri  ate  desi  gn  the  fati  gue  I  ife  of  the  damped  repai  r  can  be  i  nfi  nite,  whi  I  e  I  arge  metal  I  i  c 
repairs  operate  above  the  threshold  for  crack  growth  and  as  a  result  need  to  be  replaced 
frequently.  Over  a  ten-year  period  the  batch  of  11 F/  A-18's  suffering  from  acoustic  fatigue 
damage  in  the  aft  fuselage  would  probably  cost  $1.1  million  for  metallic  scab  repairs.  In 
comparison  the  use  of  highly  damped  repairs  could  result  in  cost  savings  of  one  order  of 
magnitude,  compared  to  metal  scab  repairs. 


10.  Conclusions 

A  generic  design  procedure  has  been  developed  using  a  number  of  closed  form  solutions. 
These  consist  of  a  simplified  dynamic  analysis  together  with  crack  bridging  theory. 
Furthermore,  the  solution  is  also  applicable  for  plates  containing  cracks  on  the  boundary, 
accounting  for  the  effect  of  residual  thermal  stresses.  The  procedure  uses  crack  growth 
equations  which  have  been  found  to  give  good  results  for  acoustic  fatigue  data.  It  has  been 
show  n  that  the  anal  ysi  s  for  crack  bri  d  gi  ng  can  pred  i  ct  the  repai  r  effi  ci  ency  w  i  th  al  I  the  pi  ate 
edges  being  either  simply  supported  or  clamped.  As  a  result,  this  generic  design  procedure 
will  enable  optimum  design  of  highly  damped  repairs  to  acoustically-induced  cracked  panels. 
It  has  been  found  that  a  highly  damped  repair  is  a  viable,  low  risk,  low  cost  repair  to  prevent 
any  crack  growth  on  the  aft  fuselage  of  the  F/  A-18.  This  repair  does  not  in  any  way  pi  ace  a 
performance  I  imitation  ontheF/  A-18.Thisrepairwould  extend  thelifeofthiscomponent  of 
theF/  A-18  beyond  6000  hours.  Furthermoreover  a  period  of  10years,thecost  savings  areone 
order  of  magnitude  improvement  than  those  for  metallic  repairs. 
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AppendixA:  Validation  Studies 

A.l.  Validation  of  un-repaired  panel 

This  section  provides  the  validation  of  the  closed  form  solution  using  both  experimental  and 
F.E.  results.  The  latter  involves  the  use  of  the  random  vibration  capability  of  theNASTRAN 
F.E.  program.  Furthermore  the  experimental  data  used  to  validate  the  closed  form  solution 
w  as  from  f I  i  ght  tri  al  s  [  31] . 


Panel  thickness  is  0.9  mm 


Figure  Al:  Dimensions  of  panel  and  dotted  outline  of  chemical  milled  area 


Figure  A  2:  F.E.  model  of  un-patched  panel 
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The  dimensions  of  pi  ate  to  be  considered  are  shown  in  Figure  A 1.  The  corresponding  F.E. 
mesh  isshown  in  FigureA2has  been  used  for  comparison  purposes.  This  geometry  is  similar 
to  a  panel  on  the  F/  A-18  aft  fuselage  where  acoustic  fatigue  cracking  has  been  found.  The 
geometry  has  been  si  mpl  if  i  ed ,  but  contai  ns  a  chemi  cal  I  y  mi  1 1  ed  regi  on .  A I  so  a  I  and  outer  area 
exists  which  is  thicker  than  the  chemically  milled  region.  This  region  is  fastened  to  the 
stiffener  with  one  or  two  lines  of  rivets.  Theouter  boundary  in  FigureA2  coincides  with  the 
rivet  lines. 


FigureA3:  Rms  microstrain  for  un- patched  panel,  versus  distance  across  the  panel 

The  results  shown  in  Figure  A3  correspond  to  the  F.E.  mesh  which  isshown  in  Figure  A2. 
These  results  show  the  comparison  of  the  closed  form  solution,  F.E.  and  flight  trial  data 
points.  The  overall  sound  pressu  re  lev  el,  0  ASP  L,  is  determined  by  matching  the  closed  form 
solution  and  the  experimental  dataatthecentreoftheplate.TheF.E.  results  determined  at  the 
calculated  OASPL  agreewell  with  the  closed  form  solution,  hence  the  experimental  data.  At 
the  ed  ge  of  the  p  I  ate  the  cl  osed  form  sol  uti  on  si  i  ghtl  y  over  esti  mates  thef  I  i  ght  tri  al  d  ata  poi  nt, 
while  the  F.E.  slightly  under  estimates  the  closed  form  solution.  Overall  the  agreement 
betw  een  th e  cl  osed  f or m,  F .  E .  an d  f I  i  g ht  d  ata  i  s  good .  A  s  a  resu  1 1  th e  I  oad  i  n g  to  ach  i  eve  th ese 
was  an  OASPL  of  159  dB.  Results  are  also  shown  in  Table  Al,  and  good  comparisons  have 
been  madefor  parameters  such  as  frequency,  lateral  displacement  and  stress  intensity  factor. 
The  value  of  the  lateral  displacement  indicates  that  a  linear  solution  is  acceptable.  The  F.E. 
model  includes  the  chemically  milled  area  together  with  the  land. 
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Table  A 1:  F.E.  validation  models  for  panel,  at  0A5PL  =159dB,  loss  factor  =  3.2% 


Model 

Frequency 

F.E. 

Frequency 

dosec 

(Hz) 

F.E. 

Lateral 

Displ. 

Closed 

Lateral 

displ 

F.E.  strain 
at  edge  of 
Panel 

microstrain 

Closed  form 
Strain 
microstrain 

Chem.  mill 

324.4 

328.8 

0.673 

0.761 

427.4 

389.5 

Skin -Hand 

374.6 

328.8 

0.692 

0.761 

340.5 

389.5 

Figure  A  4:  F.E.  mesh  density  for  chemically  milled  model  to  validate  a  closed  form  solution  for  the 
patched  model,  using  Dyad  606.  Dimensions  are  175  mm  by  176  mm. 

A  .2.  Validation  of  repaired  panel 

The  second  part  of  the  validation,  now  knowing  the  sound  pressure  level,  was  to  predict  the 
microstrain  and  compute  the  stress  intensity  for  a  repaired  plate.  The  mesh  isshown  in  Figure 
A4and  the  results  areshown  in  Figure  A  5.  Firstly  for  comparison  purposes  no  crack  exists  in 
the  reinforced  plate.  In  this  particular  case  both  the  F.E.  and  closed  form  solution  under 
estimate  the  microstrain.  Secondly  the  crack  is  introduced.  The  results  shown  in  Table  A2 
cover  the  cracked  and  un-cracked  case.  These  results  show  good  agreement  for  frequency  and 
lateral  displacement. The stressintensity factor  isover  predicted  and  henceconservativefora 
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SPL  of  159dB.  Itshould  be  noted  that  w  hi  I  ethetest  flight  value  of  the  repaired  plate  strain  at 
the  edge  of  the  plate  for  flight  1  was  87  rms  microstrain,  for  flights  2  and  4  was  83  and  72 
microstrain  respectively.  Overall  these  results  have  shown  a  good  correlation  with  testflight 
strains. 


Figure  A  5.  D  istance  across  the  repair  (A  to  B)  versus  rms  microstrain 


Table  A  2.  F.E.  validation  model  for  repaired  panel  with  crack,  2a=25.4,  0  A  SPL  =159  dB,  Dyad  606, 
load  factor  L=l. 


Temp  °C 

Freq.  F.E. 

(Hz) 

No  crack 

Freq. 

Closed 

(Hz) 

F.E.  Lateral 
Disp(mm) 
rms 

Closed 
Lateral 
Disp  (mm) 
rms 

R  peak 

F.E. 

(M  Pa  m°-5) 

Kpeak 

Closed  fonn 
(MPamO-s) 

30 

560.7 

571.5 

0.1159 

*nu  pi  atcc.dat 
*nuplatcc.db 

0.08714 

4.277 

nuplateedat 

nuplatee.db 

3.689 

50 

314.0 

314.8 

0.14218 

*nuplatdd.dat 

*nuplatdd.db 

0.1245 

3.163 

nuplataa.dat 

nuplataa.db 

2.258 

*  N  o  crack. 

F.E.  model  in  which  crack  exists  is  nuplataa.db,  nuplataa.dat,  nuplatee.db,  nuplatee.dat 


45 


DSTO-RR-0283 


TableA3:  Extension  of TableA2 


Shear 

Modulus 

(Mpa) 

Loss 

factor 

microstrain 

fmm 

fbm 

fbb 

fmb 

70 

0.115545 

72.9269 

0.1697 

5.3280E-2 

5.9360E-2 

7.8215E-3 

5 

0.394547 

68.1903 

0.18962 

3.78565E-2 

3.55668E-3 

7.8963E-2 

TableA4:  Extension  of TabieA2 


Material 
loss  factor 

**Load  factor 

L 

Program 

COD 

(mm) 

F.E.  datafile 

Crack  tip 
Size  (mm) 

0.330 

1.000 

^Damped  Repair 
16cca.nb 

1.5284E-7 

nuplatee.dat 

0.06314 

1.000 

1.000 

^Damped  Repair 
16cca.nb 

1.1298E-7 

nuplataa.dat 

0.06314 

**  load  factor  L= 1.0,  *M  al 

ihematica  file. 

(See  peak  value  in  Figure  27) 


For  the  F.E.  analysis  the  stress  intensity  factor  is  calculated  using  a  crack  opening 
displacement  formula  for  plane  stress: 


K 


rms 


EV  [2x 
4  V  / 


(Al) 


w  here  /  i  s  the  I  ength  of  the  crack-ti  p  el  ement 

1/  i  s  the  i  n-pl  ane  d i  spl  acement  at  d i  stance/  from  the  crack-ti  p 
E  is  Young's  modulus 

The  cal  cu  I  ati  on  of  the  stress  i  ntensi  ty  i  ndi  cates  the  effect  of  bou  nd  ary  cond  i  ti  ons  on  the  repai  r . 
The  difference  between  the  experimental  and  calculated  strains  is  most  likely  due  to  edge 
constraint  conditions.  The  higher  experimental  values  result  from  the  repair  extending  just 
outside  the  rivet  line.  If  the  repair  on  the  flight  trial  aircraft  were  extended  further  into  the 
adjacent  panel,  lower  edge  strainswould  havebeen  recorded.  F.E.  analysishas  shownthatit 
i  s  onl  y  necessary  to  extend  the  repai  r  across  the  I  and  regi  on  adjacent  to  the  crack. 

A  .3.  C  orrecti  on  factors 

There  is  some  disagreement  between  the  calculation  of  stress  intensity  between  F.E.  and 
closed  form  solution.  This  is  partly  due  to  the  simplification  of  the  crack  bridging  theory  in 
that  the  crack  is  normally  considered  at  the  centre  of  the  panel  rather  than  at  the  edge.  The 
F.E.  result  is  the  most  accurate.  The  stress  intensity  factor  has  been  evaluated  at  both  the 
centre  and  at  the  edge  of  the  plate  by  F.E.  and  closed  form  solution  for  crack  lengths  of 
2a=25.4  and  76.2  mm,  see  Table  A5.  For  our  case  the  value  has  been  interpolated  for 
2a=50  mm.  FI  owever  the  results  for  thecomparison  of  stress  intensity  factors  given  in  Table 
A2  do  not  justify  a  correction  factor. 
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TableA5:  Correction  factor  derived  from  F.E.  analysis,  for  stress  intensity  factor,  for  repaired  panel 
whereG=70.M  Pa 


Crack  length 
(2a) 

Centre  of  panel 

Edge  of  panel 

Ratio: 

centra  edge 

Correction 

factor 

25.4 

3.177 

cencrack.dat, 

3.2879 

nuplataa.dat 

nuplataa.db 

0.9663 

0.9663/  2=0.4831 

76.2 

3.534 

cencrack5.dat, 

cenplat83.db 

5.718 

cencrack5.dat 

cenplat86.db 

0.618 

0.618/  2=0.309 

152.4 

1.0456 

cencrack0.dat 

cenplatl58.db 

1.459 

cencrackl.dat 

cencrackl59.db 

0.716 

0.760/  2=0.358 

Interpolated  valuefor  2a=50.8  mm  gives  a  correction  factor  of  0.673. 


Figure  A  6:  D  ata  from  T  able  A  5,  crack  length  versus  stress  intensity  for  edge  crack 

The  data  from  Tabl e  A 5  for  the  stress  intensity  at  the  edge  of  the  plate  has  been  plotted  in 
FigureA6.  A  crack  length  of  2a  =50  mm  correspondsto  a  stress  intensity  of  5  M  Pa  m 05  .As 
thecrack  grows  beyond  2a=76.2  mm,  the  stress  i  ntensity  val  ueof  5  M  Pa  m 05  correspondsto  a 
crack  length  of  2a=110mm.  If  the  crack  length  of  2a=50mmissurvivable,then  beyond  a  crack 
length  of  2a=110  mm  no  crack  growth  will  occur. 


47 


DSTO-RR-0283 


Figure  A  7:  F.E.  mesh  density  for  chemically  milled  model  to  validate  a  closed  form  solution  for  the 
patched  model,  using  Dyad  606.  This  mesh  contains  a  2a=50  mm  crack. 


A  .4.  Precise  F.E.  model  for  repaired  panel 

A  large  50  mm  crack  is  shown  inthemesh  in  FigureA7.Thegeometryforthismesh  ismore 
precise  than  thatin  FigureA4,  butwith  alarger  crack.  A  traversefromA-Bfor  strain  isshown 
in  Figure  A8.  In  comparison  to  Figure  A  5  the  strain  in  the  centre  of  the  plate  has  collapsed. 
Also  the  strain  near  the  crack  is  very  low  while  the  strain  in  theplateon  the  opposite  sideof 
the  crack  is  slightly  increased.  The  stress  intensity  predicted  from  the  closed  form  solution 
agrees  w  i  th  the  crack  open  i  ng  d  i  spl  acement  resu  ItusingF.E.fortheupper  crack  ti  p .  Fi  ow  ever 
the  p red  i  cted  stress  i  ntensi  ty  for  each  crack  ti  p  are  si  i  ghtl  y  d  i  fferent.  A  change  i  n  geometry  as 
shown  inTableA6i.e.  a=173  and  a=178  mm,  does  not  significantly  change  the  result  for  the 
closed  form  solution.  Whileno  validation  for  Dyad  609  has  been  carried  out  in  Appendix  A, 
results  shown  in  Figure  27a  (Page  34)  show  very  good  agreement  in  comparison  to  F.E. 
results. 


Table  A  6:  F.E. /closed  form  solution  for  repaired  panel  with  crack,  OASPL=159  dB,  Dyad  606, 
2a=50  mm,  4  layers,  see  Figure  A  6. 


a 

(mm) 

b 

(mm) 

G 

(MPa) 

fF.E. 

(Hz) 

f  dosed 

(Hz) 

V 

LHS 

(cod) 

mm 

V 

RHS 

(cod) 

mm 

K  peak 

LHS 

F.E. 

(MPamO-s) 

K  peak 

RHS 

F.E. 

(MPam0-5) 

K  peak 

Closed 

form 

(MPa  m°-5) 

176. 

175 

70 

531. 

485. 

1.33E-7 

1.05E-7 

2.875 

3.610 

3.632 

173. 

175 

70 

508. 

- 

- 

- 

- 

3.673 

178 

175 

70 

471. 

- 

- 

- 

- 

3.598 

(Mathematica 


ile:  Damped  Repair  16.nb) 
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Figure  A  8:  D  istance  across  the  repair  (A  to  B)  versus  rms  microstrain  for  50  mm  crack 


Also  contained  i  n  Tables  A2,  A3,  A  5,  A6,  A  7  are  all  files  used  for  the  validation.  This  includes 
PATRAN  databaseses,  NASTRAN  data  files  and  Mathematica  files.  These  are  for  the  case  of 
no-repai  r  for  the  chemi  cal  I  y  mi  1 1  ed  panel ,  no  chemi  cal  mi  1 1  i  ng  and  the  repai  red  case  for  w  hi  ch 
has  only  chemical  milling.  Mathematica  users  data  shown  in  Tables  All -A23. 

Table  A  7:  Numerical  comparison  of  orthotropic  response  for  simple  test  specimen  and  contains 
stresses 


Layup 

File/  database 

(7 x  boron 

Sx 

0,90,90,0 

One  single 
element  through 
thickness 
(Homogenous) 

Valid3a 

Max=1.61E3 

Min=7.50E2 

2.89 

1.84E-2 

0,0, 0,0 

One  single 
element  through 
thickness 
(Homogenous) 

Valid  3b 

Max=1.69E3 

M  in=7.69E2 

1.88 

4.86E-2 

0,90,90,0 

Four  individual 
Elements  through 
thickness 

Valid4a 

Max=1.66E3 

Min=7.41E2 

2.89 

1.84E-2 

90,0,0,90 

One  single 
element  through 
thickness 
(Homogenous) 

Valid3a 

Max=2.97E3 

Min=1.24E2 

2.89 

1.81E-2 
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A. 5.  Simple  test  problem 

To  demonstrate  the  PATRAN  composite  I  ay  up  capability,  it  was  necessary  to  create  a  small 
test  probl  em  (Figure  A9).  This  consists  of  a  simple  beam  with  various  lay  ups  shown  in  Table 
7.  The  model  consists  of  skin,  damping  layer  and  boron  constraining  layer.  This  is  loaded  by 
an  axial  load  at  the  end  of  the  beam,  and  is  also  constrained  at  the  other.  The  results  for 
d  i  fferent  I  ay  u  ps  are  show  n  i  n  Tabl  e  A  7.  A I  so  the  cases  of  ei  ther  one  el  ement  or  fou  r  el  ements 
per  layer  are  shown. 

Results  i  n  T  able  A  7,  for  either  one  or  four  elements  per  layer,  show  a  very  smal  I  difference  in 
stresses  and  displacements.  (This  validates  the  approach  taken  for  the  acoustically  loaded 
panel  in  which  two  layers  of  boron  are  used).  Note  also  that  the  90,0,0,90  laminate  has 
i  ncreased  the  stresses.  To  repai  r  a  acousti  c  panel  w  here  the  crack  i  s  on  the  I  ongest  si  de  of  the 
panel,  the0,90,90,0  lay  up  should  beused.Theimbalancewill  be  significantly  reduced  as  the 
number  of  pi  i  es  i  ncreases  i  .e.  I  ami  nates  of  6  pi  i  es  or  greater. 


Figure  A  9:  Small  beam  with  pressure  applied  at  one  end  and  constrained  at  the  other  end 


Table  A  8:  A  ddi  ti  on  al  data  for  Table  A  7 


Layup 

4 

4 

0,90,90,0 

One  single  element 
through  thickness 

2.89 

1.84E-2 

0,0, 0,0 

Four  individual  elements 
through  thickness 

2.89 

1.84E-2 

90,0,0,90 

2.89 

1.81E-2 
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A  .6.  V  al  i dati on  of  zero  resi d ual  stresses  after  cu re 

Consider  a  circular  patch  on  a  circular  plate.  Inthiscase  the  radiusof  the  patch  is  R,  and  and 
the  plate  Ro.  In  our  case  R,  =R0  since  the  repair  overlaps  the  entire  panel  then  clamped 
boundary  conditions  are  applicable.  From  ref(52),  equ(30)  for  the  stress  under  the  repair  for 
cooling  down  after  the  cure  is  given  by: 


El 

(1-v) 


-  a{Tj  + 


(A  3.1) 


wereu  isthedisplacementattheedgeofthepatch  .  In  ourcaseu  =0.  and  duetocoolingTi  =- 
Ti  hence  we  obtain: 


crj  =  ax  T, 


(1-v) 


(A  3. 2) 


For  the  case  of  heating  the  pi  ate  alone  to  cure  temperature,  from  [32]: 


crj  =  —  0.5  aiEl 


T1  + 


(1  +  v) 

(1-v) 


TI  + 


T  -T 
1o  1 1 


2  In  (R,/R0) 


Rn 


-1 


(A  3. 3) 


This  solution  does  not  work  where  R,  =  R0  ,  however  if  R,  is  set  =  0.994i?o  then  a  good 
solution  will  be  obtained.  The  temperature  in  the  patch  is7"o  while  the  temperature  at  the 
boundary  is  7", . 

The  residual  stress  is  thus: 


<jR  =  -  0.5  axEx 


T>  + 


(1  +  v) 

(1-v) 


Tj  + 


T  -T 
1o  1 1 


2 ln(/?7  / R0) 


Ri 

Rn 


n2 


-l 


+  <2j  Tj 


Ei 

(1-v) 


(A  3.4) 


To  determine  the  residual  stress  requires  the  addition  of  equations  (A3.2  and  A3.3),  hence 
equation  (A3.4).  It  has  been  found  numerically  that  the  sum  of  stresses  cancel  out  with 
virtuality  no  residual  stress  remaining,  see  Table  A9.  Also  the  residual  stress  has  been 
computed  for  a  rangeof  poison's  ratio  and  has  been  shown  to  bevery  small  .This  exercise  has 
been  simply  to  demonstrate  that  in  this  case  residual  stresses  due  bonding  will  bezero. 


TableA9:  Residual  stresses  after  cure  from  100  °C.  Ei  =  71000  M  Pa,  ax  =23x10  6 


Ro 

(mm) 

Ro/Ri 

Ti 

°C 

To 

°C 

EquA3.1 

(MPa) 

EquA3.3 

(MPa) 

Residual  EquA3.4 
(MPa) 

88.0 

1.006 

100 

0.0 

-233.5 

233.2 

0.3 
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Table  A 10:  K/s  Dyad  609,  t  =  1.016  mm,  8  plies,  2a=50  mm.  *(D  ampedRepairl6e_609c.nb) 


Temp  °C 

G 

fmm 

fbm 

fbb 

fmb 

5.*Kmax 

L.F. 

-67 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

-60 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

-50 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

-40 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

-30 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

-20 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

-10 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

0 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

10 

300 

0.14125 

0.00622 

0.04417 

0.00625 

7.38932 

0.00867206 

20 

300 

0.14125 

0.00622 

0.04417 

0.00625 

4.6398 

0.021966 

30 

280 

0.14299 

0.00576 

0.0446 

0.00623 

3.83528 

0.0325334 

40 

200 

0.15164 

0.00341 

0.04672 

0.00607 

3.16572 

0.0512345 

50 

130 

0.16293 

0.00151 

0.04955 

0.00590 

2.36655 

0.0991672 

60 

85 

0.17415 

-0.00329 

0.05246 

0.00567 

1.97901 

0.15338 

70 

50 

0.18800 

-0.00778 

0.05628 

0.00534 

1.67625 

0.15338 

80 

25 

0.20523 

-0.01369 

0.06155 

0.00483 

1.49853 

0.335782 

90 

5 

0.23825 

-0.02547 

0.07496 

0.00349 

1.75494 

0.302044 

100 

1.7 

0.25361 

-0.03040 

0.08464 

0.00263 

2.60597 

0.147147 

110 

1.0 

0.25929 

-0.03181 

0.08952 

0.00225 

3.27513 

0.0982069 

120 

0.9 

0.26029 

-0.03014 

0.09049 

0.00297 

4.65097 

0.050208 

*Mathemati  cafile 


Table  A 11:  L  oss  factor  compon  en  ts  for  D  yad  609,  correspon  ding  T  abl  eA10,  h=1.016  mm,  2a=50  mm 


Temp  °C 

Vc 

rtx 

% 

jy 

peak 

f(Hz) 

-67 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

-60 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

-50 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

-40 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

-30 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

-20 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

-10 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

0 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

10 

0.004 

.026 

.004 

165.98 

7.0907 

813.954 

20 

0.100 

.026 

.004 

104.22 

4.6398 

814.258 

30 

0.170 

.026 

.004 

85.6162 

3.8353 

811.039 

40 

0.240 

.026 

.004 

68.7329 

3.1657 

790.063 

50 

0.400 

.026 

.004 

49.9203 

2.3666 

759.343 

60 

0.522 

.026 

.004 

40.9631 

1.9790 

721.027 

70 

0.677 

.026 

.004 

34.44 

1.6763 

663.144 

80 

0.850 

.026 

.004 

31.492 

1.4985 

573.987 

90 

1.042 

.026 

.004 

42.7501 

1.7549 

403.008 

100 

1.000 

.026 

.004 

66.9118 

2.6060 

355.326 

110 

1.000 

.026 

.004 

83.6649 

3.2751 

344.257 

120 

0.500 

.026 

.004 

117.462 

4.6510 

342.496 
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TableA  12:  Crack  bridging  data  for  Dyad  609,  t  =  0.762mm,  6  plies,  2a=50  mm. 
*(D  ampedRepairl6e_609b.nb) 


Temp  °C 

G 

fmm 

fbm 

fbb 

fmb 

Loss.Fac. 

-67 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

-60 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

-50 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

-40 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

-30 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

-20 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

-10 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

0 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

10 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0134776 

20 

300 

0.15867 

0.03457 

0.04737 

0.008913 

0.0260082 

30 

280 

0.16041 

0.03447 

0.04781 

0.008876 

0.0360016 

40 

200 

0.16890 

0.03388 

0.04999 

0.008672 

0.0540161 

50 

130 

0.17977 

0.03286 

0.05290 

0.008357 

0.100668 

60 

85 

0.19032 

0.03157 

0.05589 

0.007199 

0.154953 

70 

50 

0.20299 

0.02963 

0.05978 

0.007464 

0.240213 

80 

25 

0.21831 

0.02662 

0.06516 

0.006685 

0.354973 

90 

5 

0.24630 

0.01872 

0.07876 

0.004728 

0.368612 

100 

1.7 

0.25884 

0.01372 

0.08853 

0.003518 

0.194661 

110 

1.0 

0.26196 

0.01260 

0.09175 

0.003168 

0.133074 

120 

0.9 

0.26423 

0.01113 

0.09442 

0.0028959 

0.068694 

*Mathematicafile 


TableA  13:  Loss  factor  components,  Dyad  609,  corresponding  TableA  13,  h=0,762  mm,  2a=50  mm. 


Temp  °C 

Vc 

V\ 

M£rms 

IS 

peak 

f(Hz) 

-67 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

-60 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

-50 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

-40 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

-30 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

-20 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

-10 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

0 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

10 

0.004 

.026 

.004 

158.709 

8.3011 

735.006 

20 

0.100 

.026 

.004 

114.173 

6.2233 

735.281 

30 

0.170 

.026 

.004 

97.0556 

5.3198 

732.561 

40 

0.240 

.026 

.004 

79.9758 

4.4903 

714.646 

50 

0.400 

.026 

.004 

59.3936 

3.4111 

688.243 

60 

0.522 

.026 

.004 

49.0792 

2.8438 

654.859 

70 

0.677 

.026 

.004 

41.4192 

2.3965 

603.387 

80 

0.850 

.026 

.004 

38.0504 

2.1134 

520.961 

90 

1.042 

.026 

.004 

53.5769 

2.4318 

349.628 

100 

1.000 

.026 

.004 

85.979 

3.5892 

297.963 

110 

1.000 

.026 

.004 

108.312 

4.4346 

285.657 

120 

0.500 

.026 

.004 

151.844 

6.3139 

283.665 
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Table  A 14:  K/sDyad  609,  t  =  0.508  mm,  4  plies,  2a=50  mm.  *(D  ampedR  epai  rl6e_609a.n  b) 


Temp  °C 

G 

fmm 

fbm 

fbb 

fmb 

5*Kmax 

-67 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

-60 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

-50 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

-40 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

-30 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

-20 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

-10 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

0 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

10 

300 

0.18288 

0.078155 

0.05329 

0.013469 

10.6560 

20 

300 

0.18288 

0.078155 

0.05329 

0.013469 

8.5883 

30 

280 

0.18445 

0.078370 

0.053730 

0.013394 

7.6214 

40 

200 

0.19197 

0.079259 

0.055927 

0.012989 

6.5905 

50 

130 

0.20133 

0.079969 

0.058829 

0.012394 

5.1014 

60 

85 

0.21015 

0.080177 

0.06179 

0.011731 

4.2456 

70 

50 

0.22805 

0.079738 

0.06566 

0.010821 

3.5784 

80 

25 

0.23251 

0.078027 

0.07097 

0.009542 

3.0151 

90 

5 

0.25385 

0.069948 

0.08445 

0.006558 

3.3472 

100 

1.7 

0.26331 

0.062442 

0.09414 

0.004813 

4.8870 

110 

1.0 

0.26567 

0.059840 

0.09733 

0.004318 

6.0325 

120 

0.9 

0.26741 

0.057688 

0.09998 

0.003935 

8.5032 

*Mathemati  cafile 


Table  A 15:  Further  props.:  Dyad  609,  corresponding  TableA  16,  h=0.508  mm,  2a=50  mm. 


Temp  °C 

Vc 

lx 

% 

Loss.F 

^rms 

jy' 

peak 

f(Hz) 

-67 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

-60 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

-50 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

-40 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

-30 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

-20 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

-10 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

0 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

10 

0.004 

.026 

.004 

.0198976 

158.861 

10.2260 

647.632 

20 

0.100 

.026 

.004 

.0305986 

128.027 

8.5883 

647.859 

30 

0.170 

.026 

.004 

.0391693 

113.193 

7.6214 

645.835 

40 

0.240 

.026 

.004 

.0549778 

96.4419 

6.5905 

632.298 

50 

0.400 

.026 

.004 

.0963854 

73.8778 

5.1014 

612.196 

60 

0.522 

.026 

.004 

.146113 

61.5286 

4.2456 

586.322 

70 

0.677 

.026 

.004 

.227325 

51.8775 

3.5784 

545.341 

80 

0.850 

.026 

.004 

.346056 

47.2442 

3.0151 

476.057 

90 

1.042 

.026 

.004 

.409403 

66.6844 

3.3472 

315.022 

100 

1.000 

.026 

.004 

.233752 

108.418 

4.8970 

261.868 

110 

1.000 

.026 

.004 

.163993 

137.056 

6.0325 

248.907 

120 

0.500 

.026 

.004 

.086328 

190.856 

8.5032 

246.756 
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Table  A 16:  K/sDyad  606,  t  =  0.508  mm,  4  plies,  2a=25.4  mm. 


Temp  °C 

G 

fmm 

fbm 

fmb 

Fbb 

Loss  Factor 

-40 

300 

0.15240 

0.05753 

0.01015 

0.05042 

0.01952- 

-20 

300 

0.15240 

0.05753 

0.01015 

0.05042 

0.01953 

0 

300 

0.16583 

0.05465 

0.00841 

0.05709 

0.01953 

10 

100 

0.16583 

0.05465 

0.00841 

0.05709 

0.01998 

30 

70 

0.16907 

5.329E-2 

7.822E-2 

5.937E-2 

0.1746 

40 

17 

0.18166 

4.626E-2 

5.523E-3 

6.885E-2 

0.4195 

50 

5. 

0.18863 

0.03917 

0.00382 

0.07740 

0.3953 

60 

1.5 

0.19308 

0.03220 

0.00253 

0.08569 

0.1189 

70 

0.60 

0.19530 

0.02727 

0.00180 

0.09167 

0.0400 

80 

0.40 

0.19604 

0.02523 

0.00154 

0.09419 

0.0242 

Table  A 17:  F.E.  structural  results  for  4  ply  repair,  mesh  shown  in  Figure  A  6 


Temp 

°C 

G 

(MPa) 

COD  top 
(mm) 

COD  bot 
(mm) 

Ktop 

(MPa  m°-5) 

Kbot 

(MPa  m°-5) 

Kave 

(MPa  m°-5) 

-67 

300 

2.5679E-4 

2.0528E-4 

1.3944 

1.1231 

1.2587 

-60 

300 

2.7887E-4 

2.2294E-4 

1.5143 

1.2197 

1.3670 

-50 

300 

3.0745E-4 

2.4579E-4 

1.6694 

1.3447 

1.5071 

-40 

300 

3.3368E-4 

2.6676E-4 

1.8119 

1.4594 

1.6357 

-30 

300 

3.5799E-4 

2.8620E-4 

1.9439 

1.5658 

1.7548 

-20 

300 

3.8068E-4 

3.0434E-4 

2.0671 

1.6650 

1.8660 

-10 

300 

4.0216E-4 

3.2151E-4 

2.1837 

1.7589 

1.9713 

0 

300 

4.2248E-4 

3.3775E-4 

2.2940 

1.8478 

2.0709 

10 

300 

4.4193E-4 

3.5331E-4 

2.3997 

1.9329 

2.1663 

20 

300 

3.8943E-4 

3.1163E-4 

2.1146 

1.7049 

1.9098 

30 

280 

3.5432E-4 

2.8414E-4 

1.9239 

1.5545 

1.7392 

40 

200 

3.1802E-4 

2.5705E-4 

1.7268 

1.4063 

1.5666 

50 

130 

2.5979E-4 

2.1260E-4 

1.4107 

1.1631 

1.2869 

60 

85 

2.2609E-4 

1.8719E-4 

1.2277 

1.0241 

1.1259 

70 

50 

1.9693E-4 

1.6477E-4 

1.0693 

0.9014 

0.9854 

80 

25 

1.7646E-4 

1.4871E-4 

0.9582 

0.8136 

0.8859 

90 

5 

1.8990E-4 

1.5482E-4 

1.0312 

0.8470 

0.9391 

100 

1.7 

2.6363E-4 

2.0210E-4 

1.4315 

1.1057 

1.2686 

110 

1.0 

3.2044E-4 

2.3872E-4 

1.7400 

1.3060 

1.5230 

120 

0.9 

4.3053E-4 

3.1678E-4 

2.3378 

1.7331 

2.0355 

'Top  crack  tip  size  =  0.06714  mm,  node=20415 
*Bottom  crack  tip  size  =0.06614  mm,  node=20521 
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Table  A 18:  F.E.  structural  results  for  4  ply  repair,  mesh  shown  in  Figure  A  6 


Temp  °C 

Kp  =5*Kmax 

(MPa  m°-5) 

File 

-67 

10.8315 

rjcOLdat 

-60 

10.8315 

rjcOLdat 

-50 

10.8315 

rjcOLdat 

-40 

10.8315 

rjcOLdat 

-30 

10.8315 

rjcOLdat 

-20 

10.8315 

rjcOLdat 

-10 

10.8315 

rjcOLdat 

0 

10.8315 

rjcOLdat 

10 

10.8315 

rjc02.dat 

20 

9.5482 

rjc03.dat 

30 

8.6956 

rjc04.dat 

40 

7.8324 

rjc05.dat 

50 

6.4341 

rjc06.dat 

60 

5.6291 

rjcOLdat 

70 

4.9266 

rjc08.dat 

80 

4.4291 

rjc09.dat 

90 

4.6951 

rjcl0.dat 

100 

6.3426 

rjcll-dat 

110 

7.6146 

rjcl2.dat 

120 

10.1765 

rjcl3.dat 

rjc01.dat  contai  ns  sub  cases  of  I  oad  factors  for  PSD . 


Table  A 19:  F.E.  Temperature  results  for  4  ply  repair,  mesh  shown  in  Figure  A  6 


Temp.0 

C 

Datafiles 

(NASTRAN) 

COD  top 
(mm) 

COD  bot 
(mm) 

Ktop 

M  Pa  m0-5 

Kbot 

M  Pa  m°5 

Kther 

M  Pa  m05 

-67 

temp_87.dat 

5.5698E-3 

4.7084E-3 

30.2418 

25.7581 

28.0000 

-60 

temp_80.dat 

5.1216E-3 

4.3296E-3 

27.8083 

23.6858 

25.7471 

-50 

temp_70.dat 

4.4814E-3 

3.7884E-3 

24.3322 

20.7251 

22.5287 

-40 

temp_60.dat 

3.8412E-3 

3.2472E-3 

20.8562 

17.7643 

19.3103 

-30 

temp_50.dat 

3.2010E-3 

2.7060E-3 

17.3802 

14.8036 

16.0919 

-20 

temp_40.dat 

2.5608E-3 

2.1648E-3 

13.9041 

11.8429 

12.8735 

-10 

temp_30.dat 

1.9206E-3 

1.6236E-3 

10.4281 

8.8822 

9.6552 

0 

temp_20.dat 

1.2804E-3 

1.0824E-3 

6.9521 

5.9214 

6.4368 

10 

temp_10.dat 

6.4020E-4 

5.4120E-4 

3.4760 

2.9607 

3.2184 

20 

0 

0 

30 

temp30.dat 

-6.4165E-4 

-5.4389E-4 

-3.4841 

-2.9755 

-3.2298 

70 

temp70.dat 

-3.5460E-3 

-3.2777E-3 

-19.2545 

-17.9317 

-18.5931 

120 

templ20.dat 

-1.0475E-2 

-1.0692E-2 

-56.8806 

-58.4945 

-57.6896 
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Table  A  20:  F.E.  results  for  4  ply  repair,  mesh  shown  in  Figure  A  6 


Temp  °C 

G 

(MPa) 

Vc 

rh 

% 

Loss 

Factor. 

File 

-67 

300 

0.004 

.026 

.004 

0.019898 

rjcOLdat 

-60 

300 

0.004 

.026 

.004 

0.019898 

rjcOLdat 

-50 

300 

0.004 

.026 

.004 

0.019898 

rjcOLdat 

-40 

300 

0.004 

.026 

.004 

0.019898 

rjcOLdat 

-30 

300 

0.004 

.026 

.004 

0.019898 

rjcOLdat 

-20 

300 

0.004 

.026 

.004 

0.019898 

rjcOLdat 

-10 

300 

0.004 

.026 

.004 

0.019898 

rjcOLdat 

0 

300 

0.004 

.026 

.004 

0.019898 

rjcOLdat 

10 

300 

0.004 

.026 

.004 

0.019898 

rjc02.dat 

20 

300 

0.100 

.026 

.004 

0.030599 

rjc03.dat 

30 

280 

0.170 

.026 

.004 

0.039169 

rjc04.dat 

40 

200 

0.240 

.026 

.004 

0.054978 

rjc05.dat 

50 

130 

0.400 

.026 

.004 

0.096385 

rjc06.dat 

60 

85 

0.522 

.026 

.004 

0.146113 

rjcOLdat 

70 

50 

0.677 

.026 

.004 

0.227325 

rjc08.dat 

80 

25 

0.850 

.026 

.004 

0.346056 

rjc09.dat 

90 

5 

1.042 

.026 

.004 

0.409403 

rjcl0.dat 

100 

1.7 

1.000 

.026 

.004 

0.233752 

rjcll-dat 

110 

1.0 

1.000 

.026 

.004 

0.163993 

rjcl2.dat 

120 

0.9 

0.500 

.026 

.004 

0.086328 

rjcl3.dat 

TableA21:  Strain  for  Dyad  609,  2a=50  mm 


Temp  °C 

4  ply 

*(DR  16_e_609a.nb) 

6  ply 

*(DR  16_e_609b.nb) 

8  ply 

*(DR  16_e_609c.nb) 

-67 

152.44 

152.295 

159.271 

-60 

152.44 

152.295 

159.271 

-50 

152.44 

152.295 

159.271 

-40 

152.44 

152.295 

159.271 

-30 

152.44 

152.295 

159.271 

-20 

152.44 

152.295 

159.271 

-10 

152.44 

152.295 

159.271 

0 

152.44 

152.295 

159.271 

10 

152.44 

152.295 

159.271 

20 

128.027 

114.173 

104.22 

30 

113.193 

97.0556 

85.6162 

40 

96.4419 

79.9758 

68.7329 

50 

73.8778 

59.3936 

49.9203 

60 

61.5286 

49.0792 

40.9631 

70 

51.8775 

41.4192 

34.44 

80 

47.2442 

38.0504 

31.492 

90 

66.6844 

53.5769 

42.7501 

100 

108.418 

85.979 

66.9118 

110 

137.056 

108.312 

83.6649 

120 

190.856 

151.844 

117.462- 

*M  athemati ca  fi  I  e  i s:  Damped Repai  rl6_e_609a.nb  etc. 
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T  able  A  22:  F25_A  lum.txt,  ( 2a=50  mm  crack  length) 


Temp  °C 

fmm 

fbm 

Fmb 

fbb 

-67. 

0.10052 

0.03248 

0.008905 

0.03179 

-60 

0.10052 

0.03248 

0.008905 

0.03179 

-50 

0.10052 

0.03248 

0.008905 

0.03179 

-40 

0.10052 

0.03248 

0.008905 

0.03179 

-30 

0.10052 

0.03248 

0.008905 

0.03179 

-20 

0.10052 

0.03248 

0.008905 

0.03179 

-10 

0.10052 

0.03248 

0.008905 

0.03179 

0 

0.10052 

0.03248 

0.008905 

0.03179 

10 

0.10052 

0.03248 

0.008905 

0.03179 

20 

0.10052 

0.03248 

0.008905 

0.03179 

30 

0.10052 

0.03248 

0.008905 

0.03179 

40 

0.10052 

0.03248 

0.008905 

0.03179 

50 

0.10052 

0.03248 

0.008905 

0.03179 

60 

0.10052 

0.03248 

0.008905 

0.03179 

70 

0.10052 

0.03248 

0.008905 

0.03179 

80 

0.10052 

0.03248 

0.008905 

0.03179 

90 

0.10052 

0.03248 

0.008905 

0.03179 

100 

0.10052 

0.03248 

0.008905 

0.03179 

110 

0.10052 

0.03248 

0.008905 

0.03179 

120 

0.10052 

0.03248 

0.008905 

0.03179 

Table  A  23:  F25_Boron_c.txt ,  (cracklength  2a  =  50  mm) 


Temp  °C 

fmm 

fbm 

Fmb 

fbb 

-67. 

0.1722 

0.0732 

0.0142 

0.0397 

-60 

0.1722 

0.0732 

0.0142 

0.0397 

-50 

0.1722 

0.0732 

0.0142 

0.0397 

-40 

0.1722 

0.0732 

0.0142 

0.0397 

-30 

0.1722 

0.0732 

0.0142 

0.0397 

-20 

0.1722 

0.0732 

0.0142 

0.0397 

-10 

0.1722 

0.0732 

0.0142 

0.0397 

0 

0.1722 

0.0732 

0.0142 

0.0397 

10 

0.1722 

0.0732 

0.0142 

0.0397 

20 

0.1722 

0.0732 

0.0142 

0.0397 

30 

0.1732 

0.0735 

0.0142 

0.0399 

40 

0.1744 

0.0738 

0.0142 

0.0402 

50 

0.1755 

0.0741 

0.0142 

0.0404 

60 

0.1768 

0.0745 

0.0142 

0.0407 

70 

0.1781 

0.0748 

0.0142 

0.0409 

80 

0.1794 

0.0752 

0.0142 

0.0412 

90 

0.1808 

0.0756 

0.0142 

0.0415 

100 

0.1824 

0.0760 

0.0142 

0.0419 

110 

0.1840 

0.0764 

0.0142 

0.0422 

120 

0.1857 

0.0769 

0.0141 

0.0426 
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Appendix  B:  Costs 

B  .1.  C  osts  of  dam  ped  repai  r 

The  costs  of  this  type  of  repair  to  an  F/  A-18aftfuselagepanel  areshown  inTableBl.  The  cost 
of  a  RAAF  metal  scab  repair  are  also  shown  as  a  comparison.  Included  are  costs  of  materials 
and  labour.  Labour  isrequired  to  prepare  surfaces  of  the  damping  material  and  alsothelay- 
up  of  the  Boron,  and  adhesive,  and  placing  in  theauto-dave.  Also  the  labour  is  involved  in 
setting  up  the  repair,  prepare  the  metal  surface,  vacuum  bagging  and  bonding  and  curing  the 
patch  on  the  aircraft.  The  size  of  the  patch  for  damping  repairs  is  usual  ly  larger  than  the  panel 
i  n  order  that  the  cracks  at  the  edge  of  the  repai  r  are  covered .  I  n  thi  s  case  the  chemi  cal  mi  1 1  ed 
panel  i s  of  di mensions  176 x  176  mm.  A  factor  1.25  is  used  to  i ncreasethe patch  di  mensi on  to 
220x  220  mm.  These  costs  were  based  on  2006  costs.  Labour  to  construct  patch  was  based  on 
costs  to  prepare  4  ply  specimen  for  the  flight  trials.  A  total  cost  of  $7000.  is  slightly  I  ess  than 
the  costs  of  $10,000.  for  a  metallic  repair  used  by  the  RAAF,  which  are  replaced  every  year. 
The  damped  repair  has  an  infinite  fatigue  life. 


TableBl:  Costing  for  metal  scab  repair  and  damped  repair,  (2006  costs) 


RAAF  metal  scab  repair 

DSTO  highly  damped  repair 

ASTFITT 

20hrsfor  1  person 

Adhesive  FM300-2K 

$243 

3hrsfor  2nd  person 

Boron 

$388 

ATECH 

8hrs  for  4  people  to  R  &  1 
the  formation  light. 

Dyad  609 

$152 

ATECH 

Construction  and  cure  in 
autoclave  24hrs 

$2880 

ASURFIN 

3hrs  for  1  person  to  stri  p 
area  and  pai  nt  patch  & 
area 

Preparation  4hrs 

$480 

Application  2x8hrs 

$1920 

Clean,  paint  and  repair 
4hrs 

$480 

Total  man  hours 

62 

Total  man  hours 

48 

Assumed  rat^  hr 

$120 

Assumed  rat^  hr 

$120 

A  pproxi  mate  cost 

$7440 

Cost 

$6543 

Cost  of  heat  treatment  of 
patch  and  man  hours 

? 

- 

- 

Cap  cost  at: 

$10,000 

Cap  cost  at: 

$7,000 

Total  time 

4  days 

Total  time 

lday 

N  DI 

Difficult 

N  DI 

Better  NDI 

Lifetime 

Temporary 

Lifetime 

Infinite 
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